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FlOO MULTIVARIABLE CONTROL SYNTHESIS PROGRAM - 
EVALUATION OF A MULTIVARIABLE CONTROL USING 
A REAL-TIME ENGINE SIMULATION 

by John R. Szuch, James F. Soeder, Kurt Seldner, and David S. Cwynar 

Lewis Research Center 

SUMMARY 

Over the past several years, aircraft operational requirements have dictated the 
development of today's sophisticated turbofan engines. Classical control synthesis 
techniques have worked for the older, simpler engines. However, a need exists for a 
control synthesis procedure that can account for multiple loop interactions and can make 
use of them to optimize engine performance. One approach to solving the multivariable 
control problem is to apply the optimal control theory. The linear quadratic regulator 
is one specific area of optimal control theory that has been applied to the engine control 
problem. The FIDO Multivariable Control Synthesis program is a cooperative effort by 
the Air Force Aero- Propulsion Laboratory and the NASA Lewis Research Center aimed 
at extending earlier linear quadratic regulator work to accomplish the design, evalua- 
tion, and testing of a practical multivariable control for the FlOO turbofan engine. This 
report covers the NASA evaluation of the multivariable control logic and implementation. 
The evaluation utilized a real-time, hybrid computer simulation of the engine. Results 
of the evaluation are presented, and recommendations concernii^ future engine testii^ 
of the control are made. The results of the evaluation indicated that the ei^ine testing of 
the control should be conducted as plaimed. With minor modifications, the multivariable 
control logic and its implementation on the NASA digital computer/controller should pro- 
vide acceptable steady-state and transient performance at all flight conditions and power 
settings. It was concluded that the computer-aided, optimal control approach to design- 
ing multivariable controls does result in a practical solution to the engine control prob- 
lem. 


INTRODUCTION 


Over the past several years, aircraft operational requirements have dictated the de- 


velopment of gas turbine engines having increased performance capabilities over a 
wider operatii^ envelope. These development efforts have resulted in today’s sophisti- 
cated turbofan engines and will, undoubtedly, lead to increasingly complex ermines in the 
future. For example, future variable- cycle ei^ines may incorporate variable fan, com- 
pressor, turbine, and exhaust nozzle geometry to optimize overall engine performance 
(ref. 1). The trend toward more complex ermines has resulted in additional require- 
ments for the control system. Controls for future ermines will require measuring more 
engine variables (perhaps 10 to 20) to control engine fuel flows, and the variable geom- 
etry. The use of a closed-loop control to achieve increased control accuracy and re- 
sponse is replacing the scheduling type (open-loop) controls used on older engines. 

Classical control design (synthesis) techniques, which involve the design and eval- 
uation of single-input, single-output control loops, have worked for the older, simpler 
ei^ines. However, such techniques are cumbersome and time consiuning when applied 
to the multivariable control problem. This is due to the inherent loop interactions that 
exist but are not considered in the classical design process. A need exists for a more 
suitable control S3mthesis procedure - that is, one that can account for these loop inter- 
actions and, possibly, make use of them to optimize engine performance. 

One approach to solving the multivariable control problem is to apply the modern 
(optimal) control theory. This approach appears to be suited to the ei^ine control prob- 
lem due to the emphasis on maintaining optimum engine performance in the presence of 
a wide variety of external disturbances (i. e. , aircraft maneuvers, horsepower extrac- 
tion, etc. ). The linear quadratic regulator (LQR) is one specific area of modern control 
theory that has been successfully developed and applied to a wide variety of linear, 
multivariable control problems (ref. 2). There have also been some initial research 
and development efforts aimed at applyii^ the LQR theory to the design of controls for 
the nonlinear engine process (refs. 3 to 7). These efforts, however, have been limited 
to engine control over a narrow operating range (usually sea-level, static, standard- 
day conditions). 

The FIDO Multivariable Control Synthesis (MVCS) program, which is a cooperative 
effort by the Air Force Aero- Propulsion Laboratory (AFAPL) and the NASA Lewis 
Research Center, is aimed at extending the earlier LQR work to accomplish the design 
and testing of a ’’practical” multivariable control for a state-of-the-art turbofan engine - 
that is, one capable of operating the engine over its entire operatii^ envelope. The en- 
gine selected for the MVCS program is the Pratt & Whitney F-IOO-PW- 100(3) afterburn- 
ing turbofan. In addition to the des^ of a control for the FlOO engine, the MVCS 
program goals include the identification of advantages and disadvantages of the LQR 
method of designing ei^ne controls, the evaluation of the control design using a real- 
time, hybrid computer simulation of the engine (ref. 8), and, finally, the demonstration 
of multivariable control of the FlOO engine in a NASA altitude test facility. 
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Figure 1 illustrates the design and evaluation process used in the MVCS program. 

The program is highly dependent on cooperation and coordination between the various 
participants. The Air Force has provided the contract vehicle for the two contractors 
(i. e. , the engine manufacturer and the controls research organization) and is therefore 
responsible for monitoring the activities of the two contractors. Pratt & Whitney has 
been contracted (F33615-75-C-2048) to provide the necessary ei^ine support, while 
Systems Control Inc. has been contracted (F33615-75-C-2053) to design the control. 

Pratt & Whitney has the prime responsibility of defining the FIDO engine perform- 
ance (both steady-state and transient) by means of a nonlinear, digital computer simula- 
tion (CCD 1103-1. 0) of the engine (ref. 9). This simulation has been supplied to both 
Systems Control hic. and NASA for use in their preliminary evaluations of the control 
design. The digital engine simulation also forms the basis for the set of linear engine 
models that Pratt & Whitney has supplied to Systems Control Inc. for use in the LQR 
design procedure (ref. 9). The real-time, hybrid computer simulation of the ermine 
(ref. 8) was developed by NASA to match steady-state and transient data obtained from 
the Pratt & Whitney digital simulation. 

Pratt & Whitney has also defined the basic requirements of an FIDO control system. 
These requirements (described in CONTROL REQUIREMENTS section) form the basis 
for the control evaluation criteria used by NASA in the hybrid evaluation. By partici- 
pating in both the hybrid evaluation and subsequent engine tests, Pratt & Whitney has an 
opportunity to evaluate this new control concept using their intimate knowle(%e of the 
design problems associated with controlling gas turbine engines. 

Systems Control Inc. has the prime responsibility for developing the multivariable 
control logic (refs. 10 to 11). They must demonstrate that this Ic^ic can successfully 
control the Pratt & Whitney digital simulation of the ermine, be adequately prc^rammed 
on the NASA digital computer/controller (ref. 12), successfully control the hybrid com- 
puter simulation of the engine, and, ultimately, control the FlOO engine in a NASA alti- 
tude test facility. The controls contractor participates in both the hybrid computer 
evaluation and the engine tests to ensure that the control implementation reflects the de- 
sign intent and that the design is well understood by the other participants. 

NASA is responsible for providing the digital computer/contr oiler and the manpower 
required to (1) prc^ram the multivariable control Ic^ic on that computer, (2) develop the 
real-time, hybrid computer simulation of the FlOO engine, (3) carry out the hybrid eval- 
uation, and (4) plan and conduct the engine tests. NASA is also responsible for providii^ 
the necessary ermine -computer interface hardware. The total MVCS program involves 
a 19 -month phase I effort culminating in the hybrid evaluation and a 6 -month phase n ef- 
fort concluding with the FlOO engine tests. 

This report covers the hybrid evaluation phase of the MVCS program. Results of 
the evaluation are presented and recommendations concerning future engine testii^ of 
the control are made. Separate NASA reports will cover the implementation of the mul- 
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tivariable control logic on the NASA digital computer/controller and the results of the 
engine tests. 

Detailed discussions of the work performed by Pratt & Whitney and Systems Control 
Inc. , as part of the MVCS program, are found in references 9 to 11. 
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SYMBOLS 

cross-sectional area, m^ 
altitude, km 

fraction of compressor inlet airflow used for bleed 
discharge coefficient 

high power gain matrix, 5 x 15, appropriate units 

integral gain matrix, 5x8, appropriate units 

fan inlet guide vane angle, deg 

low power gain matrix, 5 x 15, appropriate units 

LQR gain matrix, 5x5, appropriate units 

fan discharge AP/P parameter 

low-pressure (fan) turbine inlet temperature, K 

thrust, N 

gain interpolation parameter 
aircraft Mach number 
fan speed, rpm 
compressor speed, rpm 

O 

total pressure, N/cm 
power lever angle, deg 
effective power lever ai^le, deg 

2 

total minus static pressure difference, N/cm 
(total minus static)/total pressure parameter 
compressor torque, N-cm 
fan torque, N-cm 

high-pressure turbine torque, N-cm 
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low-pressure turbine torque, N-cm 
compressor stator vane angle, deg 
Laplace operator, sec~^ 
total temperature, K 

control vector, 5 -dimensional, appropriate units 
mass flow rate, kg/sec 
fuel flow rate, kg/hr 

state vector, 5-dimensional appropriate units 
output vector, 5 -dimensional, appropriate units 
fan inlet total pressure, divided by sea level atmospheric pressure 
fan inlet total temperature divided by standard day temperature 
time constant, sec 
Subscripts: 

ADJ adjusted 

BLC compressor discharge bleed 
BLLT bleed for low-pressure turbine cooling 
BLHT bleed for high-pressure turbine cooling 
BOM hybrid with bill of material control 
com command 

DIG CCD1103-1. 0 with bill of material control 

ID fan hub 

i station number (i = 0, 2, 2.1, 2.2, 3, 4, 4.1, 5, 6, 7, 8, 13, 16) 

m measured 

max maximum 

min minimum 

MVC hybrid with multivariable control 

N net 

nom nominal, rate-limited reference -point trajectory 

OD fan tip 

RAM ram drag 


Qlt 

ROW 

S 

T 

U 

w 

X 

Y 
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ss reference-point schedule 
Superscript: 

(i)* entrance to mixing volume at station i 


ENGINE DESCRIPTION 

The engine selected for the MVCS prc^ram is the Pratt & Whitney F100-PW-100(3) 
afterburning turbofan (ref. 9). This engine is representative of current high-technolc^y 
engines and is illustrated in figure 2. The FlOO engine is a low -bypass-ratio, twin- 
spool, axial-flow turbofan. A single inlet is used for both the fan airflow and the engine 
core airflow. Airflow leaving the fan is separated into two flow streams: one stream 
passing through the engine core and the other stream passing through the annular fan 
duct. A three-stage fan is connected by a through-shaft to the two-stage, low-pressure 
turbine. A ten-stage compressor is connected by a hollow shaft to the air-cooled, two- 
stage, high-pressure turbine. The fan has variable, trailing edge, inlet guide vanes 
(CIW*s). The inlet guide vanes are positioned to improve inlet distortion tolerance and 
fan efficiency. The compressor has a variable inlet guide vane and two variable stator 
vanes (RCW's). These vanes are positioned to improve starting and to provide good 
high Mach nxmiber characteristics. Airflow bleed is extracted at the compressor exit 
and is discharged through the fan duct for startii^. Bleed is also extracted to satisfy 
installation requirements and to provide turbine coolii^. The main burner consists of 
an annular diffuser and a chamber with 16 fuel nozzles. The engine core and fan duct 
streams combine in an afterburner which consists of a diffuser and five concentric fuel 
manifolds. The afterburner discharges through a variable convergent-divergent nozzle. 
The exhaust nozzle is a balanced-beam design with an actuated divergent flap. The var- 
iable nozzle geometry provides near-optimum nozzle area, expansion ratio, and boat- 
tail drag throughout the operatii^ range. 

The current, bill of material (BOM) control system is basically hydromechanical 
with an engine -mounted, digital, electronic supervisory control. The BOM control Icgic 
sets the main burner fuel flow, fan inlet guide vanes, compressor stator vanes, after- 
burner fuel flow, and exhaust nozzle area to satisfy the pilot’s throttle command (PLA) 
and flight condition (i. e. , altitude and Mach number). In the MVCS pregram and its 
associated control design, the afterburner fuel flow is not being considered as a control 
variable due to the cost associated with providing the necessary interface hardware. 

The compressor discharge bleed is being utilized in the multivariable control, however. 

The FIOO-PW- 100(3) eigine is being made available for testiig of the multivariable 
control as part of the NASA-AFAPL Full-Scale Engine Research (FSER) pregram. 
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Among the FSER program objectives is the engine testii^ of advanced control concepts 
such as the multivariable control. 


CXDNTROL REQUIREMENTS 

One of the responsibilities of the engine manufactTmer in the MVCS program is the 
definition of the control system requirements. These requirements serve as a guide for 
the control designer and are the driving force for makii^ the control practical. These 
requirements also serve as the basis for the evaluation of the control design. Pratt & 
Whitney has provided a set of requirements for the FIDO engine control (ref. 9). These 
requirements are intended to ensure that (1) engine limits are not exceeded, (2) steady- 
state and transient specifications are met, (3) the ei^tne performance is compatible with 
aircraft requirements, and (4) the control can maintain satisfactory performance in the 
presence of component deterioration, horsepower and bleed extractions, inlet distortion, 
and afterburner ignition pulses. The following sections provide more detailed informa- 
tion on these requirements and are included as an aid in interpreting the results of the 
hybrid evaluation. 


Engine Protection 

In the current BOM control, the fan speed is limited as a function of the fan inlet 
total temperature as shown in figure 3. Normally, the fan speed would be scheduled be- 
low the limit because of performance and stability considerations. The compressor 
speed is limited to a fixed maximum value. Normally the control would maintain the 
compressor speed safely below this value since any overspeed requires at least a visual 
inspection of the high speed rotor. 

The temperatures that must be maintained within limits to protect the engine are the 
compressor discharge temperature Tg and the high-pressure-turbine-inlet tempera- 
ture T^. The compressor discharge temperature limit is based on maximum allowable 
metal temperature and the fact that turbine cooling effectiveness decreases rapidly above 
the limiting value of the temperature. The high-pressure-turbine-inlet temperature 
limit is based on maximum allowable metal temperature and turbine life considerations. 
The current control limits the low-pressure (fan) turbine inlet temperature FTIT to 
maintain the previously described temperatures within their limits. The FTIT limit is 
a function of the fan inlet temperature as shown in figure 4. In order to help meet tran- 
sient response requirements, it is allowable to exceed the FTIT limit by up to 27 kel- 
vins for a period of not more than 0. 5 second. 

To ensure structural integrity, the burner pressure P4 is also limited. The limit 
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is normally encountered only at low altitude, high Mach number conditions. 

To ensure structural stability, the fan and compressor variable geometry must be 
scheduled in accordance with the predicted blade flutter boundaries. These boimdaries 
are shown in figures 5 and 6 for the fan and compressor, respectively. Also shown are 
the current BOM van schedules. The ground rule for fan and compressor stability is 
to maintain the fan surge margin above 0. 15 and the compressor surge margin above 
0. 05. Since at some conditions this is not possible, the minimum allowable value is that 
predicted by the CCD1103-1. 0 simulation. 


Engine Performance 

Engine thrust and fuel consumption requirements are contained in the FlOO specifi- 
cation docxunent CP2903B. These specifications are classified. For the purpose of the 
multivariable control design and evaluation, the thrust and fuel consmnption goals were 
considered to be equal to the performance predicted by the CCD1103-1. 0 simvilation. 

The relation between thrust and PLA must be essentially linear and free of abrupt 
changes. 

The required performance must be maintained by the control regardless of engine 
deterioration, except where the level of deterioration is such that engine limits would be 
exceeded. Fan and compressor efficiency decrements of 1. 0 and 2. 0 percent, respec- 
tively, are typical. High-pressure turbine efficiency losses of 2. 5 percent have also 
been experienced. 

Transient thrust requirements are specified for sea-level/static, standard-day, 
uninstalled conditions. These requirements are listed in table I. With maximum horse- 
power and bleed extractions, the response times cannot exceed 125 percent of the tabu- 
lated values. Idle thrust is defined as the lowest attainable thrust at the particular flight 
condition. Intermediate thrust is the highest nonaugmented thrust, and maximum thrust 
is the highest augmented engine thrust. For thrust increments of ±13. 34 kilonewtons, 
if you start from a stabilized thrust between 25 and 45 percent of the available inter- 
mediate thrust, the time to achieve 90 percent of the thrust change should not exceed 
1. 2 seconds. For smaller thrust increments, this response criterion may be used as a 
measure of goodness. 


Aircraft Compatibility 

Engine airflow limits are set by inlet constraints. Airflow variations are restricted 
for supersonic operation to help maintain the inlet shock at a desired location. The cur- 
rent airflow limits are shown in figure 7. 
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In order to provide accessory airflow to various aircraft subsystems, the engine 
burner pressure should be maintained at a minimum level of 33. 09 newtons per 
square centimeter. At very high altitudes, the value of at intermediate thrust will 
be lower than this value. Under these conditions, no reduction in the power setting be- 
low intermediate is allowed. 

The engine should be controlled so that thrust is insensitive to variations in horse- 
power and bleed extractions up to the point where engine limits would be exceeded. For 
evaluation purposes, typical aircraft power extraction and bleed requirements are shown 
in figures 8 and 9, respectively. 


Control Sensitivity 

The engine must operate satisfactorily in the face of inlet pressure and temperature 
variations and rates of change of inlet conditions common to the operation of a highly 
maneuverable aircraft. Steady-state inlet variations of less than 1. 0 percent are typical. 
For aircraft accelerations, rates of change of 0. 103 newton per square centimeter per 
second and 1. Ill kelvins per second are representative. For aircraft decelerations, 
the rates of change may be as much as -0. 345 newton per square centimeter per second 
and -3. 89 kelvins per second. 

The control must be able to tolerate certain levels of afterburner ignition pulses. 
Figure 10 defines the maximum allowable level of afterburner ignition pressure pulses 
as a function of the fan discharge pressure Pl3- 

MULTIVARIABLE CONTROL 

A simplified block diagram of the FlOO multivariable control system is shown in 
figure 11. The control mode is basically proportional plus integral with a feed-forward 
path to provide rapid response. The proportional control action is provided by the LQR 
with the regulator gains Cp able to affect changes in all of the available control vari- 
ables U to reduce deviations in all of the state variables X relative to the specified 
reference point The state, control, and output variables used in the multivari- 

able control are as follows: 
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State variables Control variables 


Ouiput variables 


Fan speed 
Compressor speed 
Afterburner pressure 
Fuel flow 
Burner pressure 


Fuel flow 
Nozzle area 
CIW 
RCW 
Bleed flow 


Thrust 
Fan airflow 
Burner temperature 
Fan surge margin 
Compressor surge margin 


The integral control provides steady-state trimming of the engine operating point to 
satisfy performance requirements and engine limits. The integrators are controlled by 
the engine protection logic through the limit flags. The integrators limit fuel flow and 
geometry excursions to provide safe stable engine operation. 

The steady-state reference values of the state and control variables (i. e. , and 

U__) are scheduled as functions of the pilot command (PLA) and the flight condition. 

ss 

The outputs of the reference point schedules are rate limited in the transition control to 
prevent excessive deviations which could saturate the controller. The transition control 
provides a transient “model” for the regulator to follow. 

Because of the nonlinear nature of the engine process, a single set of LQR and inte- 
grator gains will not result in satisfactory closed-loop responses at all operating condi- 
tions. To provide wide-range operation, the gains are scheduled as functions of the fan 
inlet conditions and the transition value of the compressor speed as shown in figure 11. 

It should be noted that the reference point and gain scheduling problems were not ad- 
dressed in references 3 to 7. 

The current FlOO control system uses lead compensation to correct for a slow fan 
turbine inlet temperature sensor. The Systems Control Inc. design utilizes an FTIT 
“estimator” during transient operation to predict whether an overtemperature will occur. 
The estimated value of FTIT is compared with the limiting value (see fig. 3) and, if 
required, will cause the integral control to downtrim the engine before the over- 
temperature can actually occur. 

The following sections provide a more detailed description of the various elements 
that make up the FlOO multivariable control. An understanding of the control structure 
will serve as an aid in interpreting the results of the hybrid evaluation. 


Reference Point Schedules 

One of the tasks that must be accomplished in the multivariable control design is 
the calculation of steady- state, reference values for the five state variables, five con- 
trol variables, fan turbine inlet temperature, and fan discharge AP/P. These steady- 
state values are required for biasing the LQR inputs, unbiasing the LQR outputs, and 
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providing references for the trim integrators. The reference point schedules repre- 
sent a steady-state model of the engine. Where schedule outputs feed trim integrators, 
the control will force the engine to match the schedules. For this reason, the fan dis- 
charge AP/P, fan speed, fan inlet guide vane, compressor stator vane, and compres- 
sor bleed flow schedules must be sufficiently accurate to satisfy the performance re- 
quirements. Figure 12 illustrates the Systems Control Inc. approach to the reference 
point scheduling. The airflow schedule is basic to the fan match point. The scheduled 
values for fan speed Nj^, fan discharge AP/P, afterburner pressure Pg, and fuel flow 
Wp 4 are all functions of the scheduled airflow. The scheduled airflow is a function of 
the ’’effective'* power lever ai^le PLAM, the fan inlet total temperature T 2 , and the 
flight Mach number MN. The required inlet compatibility is provided by limiting the 
scheduled airflow to the allowable inlet corridor (see fig. 7). Compressor speed is also 
scheduled as a function of PLAM and T 2 . The scheduled compressor speed (corrected 
to station 2) is used to schedule both the burner pressure and the fan turbine inlet tem- 
perature. H the pilot-set PLA will result in a scheduled burner pressure below 33.09 
newtons per square centimeter, PLAM is uptrimmed to match the minimum allowable 
burner pressure. The scheduling of the fan and compressor variable vanes is accom- 
plished in the integral trim Ic^ic, which is discussed in a following section (p, 13). The 
bleed is scheduled closed. 

The reference point schedules have been simplified to reduce computer require- 
ments. This has been accomplished by using corrected parameters. To ensure that the 
control requirements are still satisfied, absolute limits are provided in the reference 
point schedules. In addition to the inlet corridor, these include maximum fan speed 
(fig. 3), maximum and minimum burner pressure, maximum compressor speed, maxi- 
mum FTIT (fig. 4), combustor blowout limits, and maximum and minimum fuel flows. 


Transition Control 

The transition control logic operates on the reference point schedule outputs to pro- 
vide a transient model for the controller to follow during changes in either the flight 
condition or the pilot command. Figure 13 illustrates the concept used to accomplish 
this. Scheduled values of rate limits are applied to rate limited, first-order lags 

driven by the outputs of the reference point schedules X and U„„. The rate limits 

ss ss 

are scheduled as functions of the fan inlet conditions ?2 and T 2 and the transition 
value of the compressor speed. The first-order lag time constants are set for each 
variable to give a smooth transition from the rate limited portion of the transient to the 
steady-state condition and to provide rapid movement of the reference point for small 
perturbation transients . The rate limits were determined from an analysis of linear en- 
gine models. The rates were established to provide the desired rates of change of 
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thrust, high-pressure turbine inlet temperature T^, and compressor surge margin. 

The rates were determined at flight conditions (altitude/Maoh number) of sea-level/ 
static, 9. 144 kilometers (km)/0. 9, 13. 72 km/0.9, and sea-level/1. 2 for low, middle, 
and high power settings. The interpolation of control gains and rates as fxmctions of 
the instantaneous flight condition and power setting is discussed in a later section 
(P. 14). 

This rate limiting approach does not provide ’’hard” protection for temperature 
and simge margins due to the dependence on linear analysis. This hard protection is 
accomplished by the engine protection logic and the trim integrators to be described 
later (p. 14). 

The transition logic also provides for an initial fuel flow and nozzle area rate 
’•kick” for gross decelerations. The magnitude of the rate kick is a function of the 
scheduled and rate-limited compressor speed reference values, the fan inlet conditions, 
and PLAM. This kick logic is also used to signal the occurrence of gross transients to 
the integrator mode control. 

The hybrid evaluation of the multivariable control covered both small perturbations 
and gross transients. The transition control logic was fully tested to determine the 
adequacy of the rates and switching logic. In addition, transition rates, designed for 
fast engine response, were also evaluated. 


LQR Control (Proportional) 

The LQR consists of a set of proportional gains which act through all of the control 
variables to reduce deviations in all of the state variables (relative to a specified refer- 
ence). The basic regulator is described by 

5 = U„om * - X„om» <« 

where 

U 5-dimensional control vector 

Unom 5-dimensional control vector (nominal trajectory) 

X 5 -dimensional state vector 

Xnom 5-dimensional state vector (nominal trajectory) 

C 5x5 LQR gain matrix 

hr 

The LQR proportional gains were calculated by solvii^ the matrix Ricatti equation 
at an operatii^ point defined by a linear, state-variable model of the engine. A per- 
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formance index (measure of goodness), which contains weightings on all state deviations 
and all control deviations, was specified. The LQR design is intended to minimize this 
index. 

An analysis by Systems Control Inc. of the approximately 30 linear engine models 
provided by Pratt & Whitney indicated that the full range of flight conditions and power 
settings could be covered by six LQR des^ns. Gain matrices were obtained at inter- 
mediate power (PLA = 83*^) for flight conditions of sea-level/static, 9. 144 km/0. 9, 

13. 72 km/0. 9 and sea level/1. 2. Idle power (PLA = 20°) gains were calculated for the 
sea-level/static and 9. 144 km/0. 9 conditions. Analyses indicated that designs at mid- 
power conditions were not required since the engine d5mamics appeared to change dis- 
tinctly depending on exhaust nozzle and low-pressure turbine choking/unchoking. 

Figvire 14 shows the form of the resulting gain matrices. Of the 25 proportional 
gain terms. Systems Control Lie. found that 9 could be eliminated (zeroed) without sig- 
nificantly affecting the closed-loop response characteristics. A later section (p. 14) 
will describe the schedulii^ or interpolation of the 6 gain matrices based on the flight 
condition and power settii^. 

The ability of the LQR to regulate the steady-state and transient performance of the 
hybrid simulation of the engine was thoroughly tested in the evaluation. Gross tran- 
sients were run to test the gain scheduling concept while small perturbations were run 
to test specific LQR designs. Afterburner ignitions and aircraft accelerations and de- 
celerations were simulated to test the regulation in the face of external disturbances. 


Trim Control (Integral) 

Integral control action is used in the multivariable control to either control the fan 
operating point or to track an engine limit. Systems Control Inc. determined that in- 
tegral control was necessary to satisfy the FlOO control requirements in the face of 
engine-to- engine variations, component deterioration, installation requirements, and 
associated reference-point scheduling errors. Integral gains Cj were calculated to 
provide trim action with a time constant of approximately 1. 0 second without affecting 
the LQR closed-loop response. The integral gains were determined at the same six 
operating points used in the LQR design. This allowed the integral gains to be included 
in the LQR matrix structure with its interpolation mechanism. Figure 15 shows the 
form of the overall gain matrix. The integral control consists of five integrators (one 
for each control variable) with each integrator driven by a combination of eight error 
terms. The basic integral control mode uses fuel flow and nozzle area to eliminate the 
fan speed deviation and the fan discharge AP/P deviation. In this mode, the engine 
limits are inactive. The CIW, RCW, and compressor bleed integrators are single- 
input integrators designed to maintain the vanes on their nominal schedules and the 
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bleed closed in steady state. Engine limits such as maximum allowable FTIT or 
minimum burner pressure will cause the appropriate limit error to be substituted for 
the fan speed error in both the fuel flow and nozzle area integrators. 

Integration is normally allowed only near steady-state conditions. The fuel flow 
integrator, however, is allowed to trace engine limits whenever they are reached. Con- 
trol saturation will cause the appropriate integrator to be clamped. Integrators are 
allowed to wind down during transients if the sign of the error term is such to reduce 
the absolute value of the integrator output. Dead-bands are applied to aU integrator in- 
puts to avoid limit cycling which might be caused by actuator hysteresis. 

The hybrid evaluation was intended to fully test the c3peration of the integral control 
imder both steady-state and transient conditions. The ability of the integral control 
modes to satisfy the steady-state performance requirements was determined. Particular 
attention was directed toward the integral tracking of engine limits. 


Gain Scheduling 

The LQR and integral control gains are scheduled in the multivariable control to 
reflect shifts in the engine dynamics as changes occur in the flight condition or engine 
power. The gain scheduling algorithm is shown in figure 16. It interpolates between 
gains that were determined at the four high -power conditions and the two low -power con- 
ditions. For a particular flight condition, a set of low- and high-power gains is com- 
puted by the gain scheduling algorithm based on fan inlet air density. The transition value 
of the scheduled compressor speed SNCTR is then used to interpolate the low- and 
high-power gains. The break between the low- and variable-power gains occurs at com- 
pressor speeds which correspond to a PLA of approximately 36°. The break between 
variable-power gains and the high-power gains occurs at a PLA of approximately 70°. 
The 36° setting corresponds to the point where the low-pressure turbine unchokes. This 
imchoking causes a significant change in the engine's dynamic behavior. 


Engine Protection 

The engine protection logic in the multivariable control provides hard limits on the 
commands to the control actuators. The engine protection algorithm is shown in fig- 
ure 17. The engine protection logic includes fuel flow limits, variable vane limits (figs. 
5 and 6), bleed limits (0. 0 to 6. 0 percent of the compressor inlet airflow), and exhaust 
nozzle area A,^ limits. The fuel flow limits include the maximum and minimum flows 
and an acceleration schedule which is a fxmction of the measured compressor speed. 

The exhaust nozzle area limits are functions of the PLAM. The maximvim and min im u m 
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allowable area commands converge at idle power so as to prevent limit cycling in this 
sensitive operatii^ regime. This feature was included as a result of early results from 
the hybrid evaluation. Flags are set when any unlimited commands exceed the specified 
limits or when control saturation is detected. The flags are used to clamp the appro- 
priate trim int^rators. 


FTIT Estimator 

The current FIDO BOM control system uses lead compensation to predict an im- 
pending FTIT overshoot. This is necessary because of slow FTIT sensor dynamics. 
High downtrim gains are used to avoid excessive temperature overshoots, hi the multi- 
variable control, however, another method of estimating FTIT is used. As shown in 
figure 18, the FTIT "estimator” uses a combination of the sensed FTIT (FTIT^), the 
steady-state reference value of FTIT (FTIT_„), the transition value of fuel flow 
(Wf 4 nom^’ commanded fuel flow (iVp 4 to predict the final value of 

FTIT during a transient. The predicted value is compared with the maximum allowable 
FTIT (fig. 4). If an overtemperature is predicted, the fuel flow integrator downtrims 
the fuel flow before an overtemperature can occur. 

A detailed description of the FlOO multivariable control logic can be found in ref- 
erence 11 . 


CONTROL IMPLEMENTATION 

The FlOO multivariable control logic was implemented on the SEL810B digital com- 
puter (ref. 12) shown in figures 19 and 20. The salient features of that computer are 
listed in table n. Althov^h this computer is not flight-qualified hardware, its memory, 
speed, and word size are believed to be representative of the computers that will be 
used to control engines in the 1980’s. All of the SEL810B programming was performed 
in assembly language in order to reduce the core requirements and computation time. 
The total software package, including the control a^orithms, matrix data, reference 
point schedule data, and general-purpose input/output subroutines consumed approxi- 
mately 12 000 words of core storage. The multivariable control algorithm consumed 
approximately 7000 words of core storage. A complete breakdown of the memory allo- 
cation is given in table IH. The core requirements listed in table III are not a direct 
measure of the control complexity, however. While the schedules and matrix data 
would consume the same amount of core on any 16 -bit computer, the control algorithm 
requirements are highly dependent on the computer architecture, instruction set, and 
objectives of the programmer. In the MVCS program, a determining factor was the 
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stated goal of achieving a 10-millisecond control update time, as compared with the 
14-millisecond update time for the FlOO electronic supervisory control and the 20- 
millisecond update time for the IPCS (ref. 13) implementation of the TF30 BOM con- 
trol logic. To achieve the desired 10 -millisecond update time, it was necessary to 
write the entire program using inline code. This had the effect of eliminating the over- 
head associated with multilevel subroutine calls and indexii^. While speeding up the 
computations, this required additional core storage. 

In addition to performing the multivariable control calculations, the SEL8 lOB was 
used to simulate the control actuators and for collecting transient data. All of these 
tasks were performed in the 10-millisecond update interval. 


EVALUATION OF MULTIVARIABLE CONTROL 

The principle objective of the MVCS program is to demonstrate the multivariable 
control of the FlOO engine in a NASA altitude test chamber. To that end, the multivar- 
iable control logic was developed by Systems Control Inc. and implemented on the NASA 
digital computer/controller. The implementation was described in the precedii^ sec- 
tion. Prior to the altitude tests, however, it is necessary to evaluate the control logic 
and its implementation to ensure safe, stable ei^ine operation throughout the FlOO 
flight envelope. 

This evaluation was accomplished usii^ a real-time, hybrid computer simulation 
of the FlOO engine (ref. 8). The hybrid computer simulation provided a convenient and 
economical means of evaluating both the steady-state and transient performance of the 
multivariable control and its implementation. The simulation did not, however, provide 
sufficient precision and repeatability for determining the control’s sensitivity to instal- 
lation and engine component deterioration. For those evaluation items, the CCD1103-1.0 
simulation was run by Systems Control Inc. using compatible FORTRAN coding of the 
multivariable control logic. The results of that study will be presented in their final 
report. 

Figure 21 illustrates, schematically, how the various computers were configured 
for the hybrid evaluation. An Electronic Associates Inc. model 690 hybrid computer 
provided the real-time engine simulation. The analog signals, representing the engine 
variables to be sensed, were trmiked to a TR48 analog computer which provided real- 
time sensor simulations. The simulated sensor outputs were digitized and processed 
in the SEL810B digital computer. The SEL810B computer was used to compute the com- 
mands to the control actuators and to simulate, digitally, the actuator dynamics. The 
digital actuator outputs were converted to analog signals for input to the engine simula- 
tion. The following sections provide a more detailed description of the engine, sensor, 
and actuator simulations. 
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Real-Time Engine Simulation 


A real-time, hybrid computer simulation of the FlOO-PW- 100(3) turbofan engine 
(ref. 8) was developed by NASA to support controls research programs involving that 
engine. The simulation has both wide-rai^e steady-state and transient computing capa- 
bilities. The mathematical model describing the performance of the FlOO engine was 
patterned after the CCD1103-1. 0 digital simulation (r^f. 9). Wide-range, overall per- 
formance maps for the engine's rotating components (i. e. , fan, compressor, turbines) 
were utilized so as to provide the wide-range steady-state accuracy. Factors such as 
fluid momentmn, mass and energy storage, and rotor inertias were included to provide 
transient capability. 

Figure 22 contains a computational flow diagram of the real-time simulation. It 
can be seen that the mathematical model consists of a number of individual elements and 
their related volumes, each of which requires a munber of input variables and generates 
one or more output variables. However, xmlike a digital simulation, all calculations are 
considered to be performed in parallel. 

The equations describing the FlOO mathematical model were implemented on the 
Lewis Research Center’s hybrid computing system. The Electronic Associates Inc. 
model 690 hybrid computer consists of a model 640 digital computer, model 680 analog 
computer, model 681 analog computer, and an interface unit which allows communica- 
tion between the individual computers. A photograph of the hybrid computer system is 
shown in figure 23. The digital portion of the hybrid computer was used primarily to 
perform the bivariate function generation associated with modeling the performance of 
the engine components. In addition to this function generation, the digital computer was 
also used for computing the fan and compressor surge margins and the engine net thrust. 
The digital portion of the simulation consumed 12 440 words of core storage and had an 
update time of approximately 7. 5 milliseconds. The digital portion of the simulation 
was structured to provide accurate, real-time dynamics with this update time (ref. 8). 
The remaining calculations were performed on the analog computers. The analog com- 
puters provided continuous integration with respect to time, multiplication, division, 
imivariate function generation, etc. The two analog computers were fully utilized. For 
example, the full complement of 54 multipliers and a total of 189 potentiometers were 
required. The use of peripheral equipment such as X-Y plotters and strip-chart re- 
corders allowed continuous monitoring of computed variables. Setup and checkout of the 
analog portion of the simulation was accomplished usii^ the digital computer and 
general-purpose software. 


Sensor and Actuator Simulations 

To demonstrate multivariable control of the real-time engine simulation, it was 
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necessary to simulate, in real time, the dynamic characteristics of the engine- control 
interfaces. This included the sensors and the control actuators. Since the hybrid 
evaluation was aimed at verif 3 dng the multivariable logic and its implementation prior 
to engine testing of the control, it was concluded that the interface simulations should 
reflect the hardware to be used in the altitude tests. In most cases, the dynamic char- 
acteristics of the research actuators and sensors differ from those of the BOM control 
interface hardware. 

Table IV lists the simulated sensors and their dynamic characteristics. With the 
exception of the FTIT and the fan discharge temperature Tj^g sensors, the sensors 
were represented by first-order lags. The time constants for the rotor speed sensors 
were based on the CCD1103-1. 0 simulation values. The time constants for the P2, P4> 
and Pg sensors were based on the use of strain-gage-type transducers with fairly short 
sensing line lengths. Because of the uncertainty in estimating the dynamic characteris- 
tics of the research temperature instrumentation, the BOM sensor models were used 
in the hybrid evaluation. 

In the altitude tests of the control, measurements of the fan discharge total pres- 
sures and total-minus -static pressure differences (AP) will be used to compute the Sys- 
tems Control Inc. fan discharge AP/P parameter DPE. The following equation is to 
be used: 


DPE = 


(AP)i3 + (AP)2.2 

^13 ^ 2 . 2 


( 2 ) 


Since similar transducers are to be used for the four separate measurements, it was 
assumed that they would have similar dynamic characteristics. Since the ratio would 
effectively cancel their dynamic effects, no sensor djmamics were considered in the 
computation of DPE. In the hybrid computer simulation of the engine, static pressures 
are not calculated. However, the individual AP/P parameters are computed as fvmc- 
tions of the corresponding corrected airflows. In the hybrid evaluation, the followii^ 
method of computing DPE was used: 

DPE - <^P/P>13Pl3^<^P''P)2.2P2.2 (3 

Pl3+*'2.2 

The sensor dynamics were simulated in real time using an Electronic Associates model 
TR48 desk-top analog computer which is shown in figure 24. 

The actuators associated with the five control variables had to be simulated in real 
time also. Because of equipment limitations, the SEL810B was used for the actuator 
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simulations. 

The exhaust nozzle and CIW actuators were designed by NASA to support fan 
flutter studies that were conducted as part of the FSER program. The exhaust nozzle 
interface consists of an electrohydraulic servo which positions, by means of a cable 
linkage, an air motor control valve on the engine -mounted actuation system. This is 
shown schematically in figure 25. Experimental data indicated that the response of the 
electrohydraulic servo was fast compared with the actuation system slew rate. There- 
fore, the dynamics of the servo were neglected. The remainder of the actuation system 
was simulated in the SEL810B computer. 

The CIVV actuator consists of an electrohydraulic servo which drives dual engine- 
mounted piston-in-cylinder actuators. The CIW system is shown schematically in fig- 
ure 26. For the multivariable control tests, the servovalve has been sized to provide 
the desired CIW slew rate. Prior to the hybrid evaluation, the servoamplifier gain 
was adjusted to provide satisfactory simulation response characteristics. 

NASA must also provide the fuel flow, RCW, and bleed flow interface hardware 
for the multivariable control tests. The Bendix Energy Controls Division was contracted 
(NAS3-20245) to provide the fuel flow and RCW interface hardware. The fuel flow 
actuator is shown schematically in figure 27. This system consists of a basic FlOO 
metering valve positioned by a fuel-actuated power piston. This piston is driven by an 
electrically controlled servovalve. A resolver is attached to the metering valve shaft 
to provide metering valve position feedback to the servoamplifier. The pressure drop 
across the metering valve is sensed by the FlOO fuel pump controller and is maintained 
at 41. 37 newtons per square centimeter. The pump controller dynamics are represented 
by a first-order lag with a time constant of 0. 1 second. This lag was simulated on the 
analog computer while the remainder of the fuel actuation system was simulated in the 
SEL810B computer. The nonlinear resolver characteristic was included in the simula- 
tion. 

The RCW actuator is shown schematically in figure 28. This system consists of 
dual, fuel-powered pistons driven by an electrically controlled servovalve. The servo- 
valve is designed to drive both the tandem cylinder and the RCW slave cylinder. A 
resolver is attached to the tandem cylinder to provide position feedback to the servo- 
amplifier. The RCWT resolver characteristic is essentially linear and was assumed to 
be linear in the actuator simulation. Based on the CCD1103-1. 0 simulation, RCW hys- 
teresis of ±1. 0° was included in the SEL810B simulation. 

Prior to the hybrid evaluation no information was available on the bleed flow valve 
characteristics. For the evaluation the bleed system was assumed to be fast acting with 
the bleed flow computed from 


w 


BLC 
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where BLC is the fraction of the compressor inlet airflow used for bleed in the multi- 
variable control. In the altitude tests a dynamic measurement of the compressor inlet 
airflow will not be available. An analysis of CCD1103-1. 0 data has shown that the re- 
quired bleed valve opening can be related to the BLC command by 

A 0.01774 BLC _2 ,, 

^blc = ~7 “ 

'^d, BLC 


Procedure 

Table V lists the flight conditions and types of tests covered in the evaluation. The 
flight conditions were selected to cover the full range of FlOO operation. They included 
conditions to be run in the altitude facility and, additionally, low altitude and high inlet 
temperature conditions that cannot be run in the facility. The range of power settings 
at each flight condition was limited so as to satisfy the inlet airflow requirements (fig. 7). 

The EAI640 and SEL810B digital computers provided digital listings of steady-state 
engine and control data, respectively. A total of 192 ermine and control variables were 
listed at each of the 56 evaluated operating points. 

The usual peripheral analog equipment was used to monitor transient data during 
the evaluation. Because of the vast amount of available data, however, a more conven- 
ient means of recording and processing transient data was desired. Therefore, the 
SEL810B digital computer was also utilized as a transient data sampler and storage de- 
vice as shown in figure 21. A total of 72 variables were sampled and stored for each of 
the 77 transient runs. These variables included engine variables, sensed engine vari- 
ables, internal control variables, actuator commands, and actuator outputs. Two hun- 
dred samples per variable per run were stored with the sampling rate adjusted to match 
the duration of the transient. For example, each of the 72 variables was sampled every 
100 milliseconds during a 20-second run. After the completion of a particular transient 
run, the stored data were transferred to a disk storage device for later processing. 

Prior to taking data at a particular flight condition, it was necessary to "trim” the 
engine simulation to match CCD1103-1. 0 data (ref. 8). The fan and compressor dis- 
charge specific heats and the main combustor efficiency were adjusted (if necessary) to 
match the baseline values of fan speed, compressor speed, and FTIT. During this 
trimming process, the control inputs to the hybrid simulation were fixed at their 
CCD1103-1. 0 intermediate thrust values. No trimming of the simulation was done at 
part-power conditions. 

In order to avoid transients when switching to the multivariable control, the follow- 
ii^ procedure was used. After the trimming of the engine simulation, the multivariable 
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control was put in a "TRACK" mode whereby the control’s trim integrators were con- 
tinuously updated to match the open-loop control inputs to the engine. The integrators 
were then clamped and the loops between the control implementation and the engine simu- 
lation were closed. The trim integrators were then released without an accompanying 
transient. This procedure proved successful and will be used in the altitude tests when 
switching from the BOM mechanical controls to the multivariable control. 

During the initial phase of the evaluation, a number of anomalies were discovered 
which required modifications to the control logic and its implementation. With the as- 
sistance of the control designer, causes were identified and corrective action was taken. 
Table VI lists the symptoms, causes, and required control modifications. In general, 
the observed behavior was associated with the nonideal hybrid simulation environment. 
The control logic and implementation were changed to accommodate noisy signals and 
observed differences between the hybrid and CCD1103-1. 0 simulations. 


RESULTS AND DISCUSSION 
Steady-State, Uninstalled Performance 

The steady-state, uninstalled performance of the FlOO engine simulation under mul- 
tivariable control was evaluated at those flight conditions listed in table V. Figures 29 
to 35 show the results of that evaluation. 

For each flight condition, the observed values of fan speed, compressor speed, 
burner pressure, afterburner pressure, fan turbine inlet temperature, and fan discharge 
AP/P are plotted against the corresponding scheduled values over the allowable range 
of power settings. Depending on the integral trim mode, a matching of the reference 
point schedules may be forced by the trim integrators. A match of the reference point 
schedules for uncontrolled variables, however, indicates that the schedules represent an 
accurate, steady-state model of the simulated engine. On the other hand, large devia- 
tions from the reference point schedules can result in the saturation of the regulator 
and/or the integral controller and produce degraded steady-state performance. 

For each flight condition, the observed values of net thrust, fan surge margin and 
compressor surge margin are plotted against PLA. However, because of possible hy- 
brid simulation errors caused by simplified thrust and surge margin calculations, a 
direct evalviation of those aspects of the control performance was not possible. An at- 
tempt was made to estimate those errors at each operatii^ point and to adjust the simu- 
lation data accordingly. Hybrid simulation values of fan surge margin, compressor 
surge margin, and net thrust, obtained with an SEL8 lOB implementation of the BOM 
control logic, were compared with the corresponding CCD1103-1. 0 data. Any observed 
differences were attributed to the hybrid simulation and were applied to the multivariable 
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control results. The ’’adjusted” values of the surge margins and net thrust were thus 
computed from 

^ADJ " ^^DIG ■ ^BOM^ ^MVC 

The resultant data were then judged relative to the previously described control require- 
ments. It should be remembered, however, that the primary objective of the evaluation 
was to verify the control implementation and to identify control problems prior to the 
engine tests. A quantitative evaluation of the steady-state, closed-loop performance was 
a secondary objective that could be accomplished using the CCD1103-1. 0 simulation. 

The hybrid simulation is sufficiently accurate, nonetheless, so as to provide a realistic 
test of the multivariable control logic as well as the implementation. 

Figures 29 and 30 show the steady-state results for the sea-level/static and 
3. 048 km/0. 6 flight conditions, respectively. The 3. 048 km/0. 6 condition will serve 
as a substitute for the sea-level/static condition in future ei^ine tests since simulatii^ 
altitudes lower than about 3 kilometers in the NASA altitude test facility is not possible. 
At these two flight conditions the control action was quite similar. At intermediate 
power (PLA = 83°), the trim integrators for fuel flow and nozzle area maintained FTIT 
at its limit and the fan discharge AP/P parameter at its scheduled value. This is 
shown in figures 29(e) and (f) for the sea-level/static condition. The observed agree- 
ment between the other engine variables and their schedules indicated that the reference 
point schedules did match the engine simulation quite well at these conditions. At part- 
power settings, the trim integrators maintained the fan speed and the fan discharge 
AP/P parameter at their schedules values. This is shown in figures 29(a) and (f) for 
the sea- level/static condition. However, at PLA settings of 30° and below, the nozzle 
area was constrained by the multivariable control to be 0. 274 square meter to satisfy 
the idle thrust requirement and to prevent limit cycling in the low power portion of the 
operating line. Therefore, only integral control of the fuel flow was available to main- 
tain the fan speed at its scheduled value. The area constraint was active at each of the 
evaluated subsonic conditions below the 13. 72-kilometer altitude. Figures 29 and 30(g) 
show that the multivariable control did satisfy the intermediate thrust requirement at 
the three flight conditions. The idle thrust results varied, however, with a match of 
the CCD1103-1. 0 idle thrust at the sea-level/static condition, but a higher thrust was 
observed at the 3. 048 km/0. 6 condition. 

The control requirements specified minimum allowable fan and compressor surge 
margins of 0. 15 and 0. 05, respectively. These requirements are relaxed at those 
operating points where the CCD1103-1. 0 simulation predicts lower values. In these 
two cases, however, the multivariable control produced slightly lower surge margins 
than the CCDl 103-1. 0 simulation. At the sea-level/static condition, a minimum fan 
surge margin of 0. 074 was observed as shown in figure 29(h). This was slightly lower 
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than the CCDl 103-1. 0 value of 0. 079. Also, at the 3. 048 km/0. 6 condition, the multi- 
variable control resulted in a minimxun fan surge margin of 0. 134 as compared with the 
baseline value of 0. 144. This is shown in figure 30(h). 

Two different aspects of the multivariable control were introduced and verified at 
the 9. 144 km/0. 9 and 13. 72 km/0. 9 conditions. Those results are shown in figures 
31 and 32, respectively. For these conditions, the trim integrators for fuel flow and 
nozzle area maintained both the fan speed and fan discharge AP/P parameter at their 
scheduled values at the intermediate thrust setting. At this setting, the FTIT was 
lower than the limit. This is shown in figures 31(a), (e), and (f) for the 9. 144 km/0. 9 
condition. At idle power the trim integrators switched from fan speed control to burner 
pressure control so as to maintain a minimum burner pressure of 33. 09 newtons per 
square centimeter. This is shown in figures 31(a) and (c) for the 9. 144 km/0. 9 condi- 
tion. Figures 32(a) and (c) show that the minimum burner pressure limit was active at 
a higher power setting due to the higher altitude. PLAM was uptrimmed to 22. 7° at the 
9, 144 km/0. 9 condition and to 38. 1° at the 13. 72 km/0. 9 condition. Figure 31(f) shows 
that the nozzle area constraint was also active at low power settings at the 9. 144 km/0. 9 
condition. 

Figures 31 and 32(g) show that the multivariable control again satisfied the inter- 
mediate thrust requirement. In both cases, however, the idle thrust was lower than that 
predicted by the CCD1103-1. 0 simulation. This was caused by the BOM control limit- 
ii^ the burner pressure to a higher value than 33. 09 newtons per square centimeter. 

The surge margin requirements were satisfied at both flight conditions as shown in fig- 
ures 31(h) and (i) and 32(h) and (i). At the 13. 72 km/0. 9 condition, the multivariable 
control matched the CCDl 103-1. 0 minimum fan surge margin of 0. 142. 

Figures 33 and 34 show the results from the evaluation of the multivariable control 
at two low altitude, high Mach number conditions. The sea- level/1. 2 condition produced 
the highest engine pressures and flows. The 3. 048 km/1. 2 condition will be a substitute 
for the sea-level/1. 2 condition in the altitude tests with the engine. F^ures 33(a) and (f) 
and 34(a) and (f) show that the integral trims for fuel flow and nozzle area maintained the 
fan speed and fan discharge AP/P at their scheduled values over the entire operating 
line. PLA settings below 50° were not run since the control criteria do not permit 
lower settings due to inlet airflow requirements at these supersonic flight conditions. 

The results of the evaluation showed that the multivariable control did maintain the in- 
termediate operating point below the maximum burner pressure and maximum FTIT at 
at the sea-level/1. 2 condition. The resultant thrust exceeded the CCD1103-1. 0 value at 
intermediate power but was considerably lower at part-power settings as shown in fig- 
ure 33(g). Figure 34(g) shows that the intermediate thrust at the 3. 048 km/1. 2 condition 
was also low. The surge margin requirements were satisfied at both flight conditions. 

At the sea-level/1. 2 condition, the intermediate value of the compressor speed, 
shown in figure 33(b), exceeded the maximum allowable value. While this was attributed 
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to a simulation anomaly at this operating condition, it does point up the need for some 
form of hard limit on compressor speed in addition to the limit provided in the refer- 
ence point schedule. The steady-state portion of the evaluation was performed with the 
acceleration schedule of fuel flow removed from the engine protection logic. Since this 
schedule has a high and low compressor speed cutback feature, it is felt that this sched- 
ule will provide sufficient protection against a compressor overspeed during engine test- 
ing. 

Figure 35 shows the steady-state results for the evaluation points where part-power 
operation was not permitted. At these conditions, only intermediate power settings 
were run. For all but the 12. 19 km/2. 2 condition, the integral trims maintained FTIT 
at its limit. At the 12. 18 km/2. 2 condition, the fan speed was maintained at its sched- 
uled value. The fan discharge AP/P parameter was maintained at its scheduled value 
at all of the flight conditions. 

At three of the flight conditions, however, the authority limits on the fuel flow and 
nozzle area trim integrators had to be relaxed to match the schedtiles. At the 
6. 096 km/1. 8 condition, the fuel flow uptrim limit was increased to permit 878. 6 kilo- 
grams per hour of uptrim to match the FTIT limit at that condition. At the 12. 19 km/ 
2. 2 condition, the fuel flow uptrim limit was increased to permit 838. 2 kilograms per 
hour of uptrim to match the fan speed schedule. Finally, at the 19. 81 km/2. 5 condition, 
the nozzle area downtrim limit was increased to 0. 073 square meter to match the fan 
discharge AP/P schedule. 

At the 19. 81 km/0. 9 and 22. 86 km/1. 8 conditions, conflicting control requirements 
existed. That is, the burner pressure was lower than the minimum allowable value 
while the FTIT was at its limit. The integral control logic places FTIT control on a 
higher priority. Therefore, the burner pressure was lower than 33.09 newtons per 
square centimeter at these conditions as shown in figure 35(c). 

A problem with the FTIT estimator logic was discovered at the 19. 81 km/0. 9 con- 
dition. Figvire 35(e) shows a large discrepancy (140°) between the maximum and sched- 
uled values of FTIT at that condition. It was observed that changes in the fan inlet 
conditions due to changes in the flight condition or noise could produce sufficient changes 
in the scheduled value of FTIT so as to reinitialize the estimator. The estimator was 
designed to reset to the scheduled value when a transient was indicated by a sudden 
change in the schedule output. At the 19. 81 km/0. 9 condition, the resettiig of the esti- 
mator to the low value and the low burner pressure resulted in a fuel flow uptrim and a 
corresponding increase in the FTIT. An overtemperature could not be prevented by the 
existing estimator logic. Modifications to the FTIT reference point schedule to better 
match the engine at this flight condition could not conveniently be made. However, the 
estimator logic was modified to prevent the reset action when low burner pressure condi- 
tions existed. 

Figure 35(g) compares the actual net thrust obtained with the multivariable control 
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and the adjusted net thrust. While all of the flight conditions resulted in thrust in excess 
of the CCD1103-1. 0 values, only the 19. 81 km/0. 9 condition resulted in adjusted thrust 
meeting the thrust requirement. This could be attributed to the fact that no supersonic 
conditions were considered in the establishment of the reference point schedules. Fig- 
ure 35(h) shows that the surge margin requirements were satisfied at all of the flight 
conditions. The 19. 81 km/0. 9 condition produced a minimum fan surge margin of 
0. 137 which exceeded the CCD1103-1. 0 value of 0. 126. 

In general, the results of the steady-state evaluation of the multivariable control 
indicated that the control logic and its implementation did satisfy most of the control 
requirements. Engine limits were maintained throughout the flight envelope. The un- 
certainties associated with the simulation calculation of net thrust and surge margins 
made it difficult to make a quantitative assessment of those aspects of the control per- 
formance. It appears that those instances where thrust and surge margin requirements 
were not satisfied do not warrant changes to the control design. 

The reference point schedules appeared to accurately represent the engine simula- 
tion at subsonic flight conditions. In spite of differences between the hybrid and 
CCD1103-1. 0 simulations, the control provided good steady-state regulation of the 
operating point; thus, the control should be able to accommodate expected engine-to- 
engine variations. 

Although thrust requirements were specified only at low altitude, low Mach number 
conditions, it appears that some modifications to the reference point scheduled might be 
required to match CCD1103-1. 0 values of thrust at supersonic flight conditions. Those 
modifications would probably eliminate the trim saturations observed at three of the 
supersonic flight conditions. 


Power Lever Transients 

The transient response characteristics of the FlOO real-time simulation under mul- 
tivariable control were investigated at those flight conditions where part-power settings 
were permitted. Those conditions are listed in table V. The transients included small 
and large amplitude PLA snaps and cyclic movements of the power lever. For all but 
the cyclic movements, the PLA was ramped at a constant rate of ±126° per second to 
match the altitude facility rate limit. As described in preceding sections, the engine 
and control variables were sampled by the SEL810B computer and the resultir^ data 
were stored on disks for later processing. 

Figures 36 to 45 contain plots of selected transient data. Figures 36 and 37 show 
the responses of engine and internal control variables, respectively, to a power lever 
snap from 30 percent thrust (PLA = 35°) to intermediate thrust (PLA = 83°) at the 
3. 048 km/0. 9 condition. The correspondii^ deceleration is also shown. Figures 36(b) 
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to (f) show the responses of the five control variables. As shown in figure 36(b), the 
fuel flow was not constrained by the acceleration schedule for this transient. As the 
intermediate thrust was approached, a dip in the fuel flow was observed. This was fol- 
lowed by a slow increase to its final value. This behavior was attributed to the switch 
from the untrimmed, regulator control mode to the integral trim mode as the steady- 
state condition was approached. Figure 36(f) indicates that the compressor discharge 
bleed was utilized during the deceleration. A peak outflow bleed of 1. 09 kilograms per 
second was observed. This represented about 1. 7 percent of the compressor inlet air- 
flow. Figures 36(g) to (o) show the resultant responses of selected engine variables. 

As shown in figures 36(g) and (h), the fan and compressor rotor speed responses ex- 
hibited no overshoot during the acceleration and only a slight compressor overshoot dur- 
ing the deceleration. Figures 36(i) to (1) show that the burner pressure, FTIT, after- 
burner pressure, and net thrust exhibited the same characteristics as the fuel flow 
response. The responses of the FTIT sensor and estimator are also shown in fig- 
ure 36(j). At the 30 percent thrust condition, the estimator matched the sensor output. 
During the acceleration, the estimator led the actual temperature response and exceeded 
the FTIT limit by 40 kelvins, 0. 5 second after the initiation of the transient. This re- 
sulted in downtrimming the fuel flow and the subsequent limiting of the FTIT peak 
value. Figvure 36(^) shows that no thrust overshoots were observed during the accelera- 
tion and deceleration. The response time (time to achieve 90 percent of the thrust in- 
crease) of 2. 1 seconds was acceptable. Figures 36(m) and (n) show the responses of 
the fan and compressor surge margins, respectively. No loss of fan surge margin was 
observed during the acceleration. A minimum value of 0. 15 was observed for the com- 
pressor surge margin during the deceleration. This did not violate the 0. 05 minimum 
requirement. Figure 36(o) shows the response of the fan discharge AP/P parameter. 

In steady-state, this parameter was maintained at its scheduled value by the trim inte- 
grators. During the acceleration, a maximum deviation of 0. 013 was observed while the 
trim integrators were clamped. A smaller deviation of 0. 007 was observed during the 
deceleration. 

While the data presented in figure 36 provide a measure of the multivariable control 
performance, it was desirable to monitor internal control variables to verify the control 
Ic^ic. Figure 37 shows plots of LQR outputs, integral trims, and state deviations for 
the same 3. 048 km/0. 9 transient shown in figure 36. Figures 37(b) and (c) show the 
responses of the LQR and integral control outputs of fuel flow, respectively. Durii^ the 
initial portion of the acceleration, the negative fan speed error (relative to the sched- 
uled value) would have produced a fuel flow uptrim by means of the integrator. The 
multivariable control logic, however, clamped the integrator during this portion of the 
transient as shown in figure 37(c). The LQR did, however, uptrim the fuel flow to re- 
duce the state deviations as shown in figure 37(b). About 0. 5 second after the initiation 
of the transient, the high FTIT estimator output caused the fuel flow integrator to be 
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released to produce a fuel flow downtrim. The decreasing state deviations and the inte- 
gral downtrim of fuel flow resulted in a reduction in the LQR output. At about the 1. 5- 
second mark, the falling FTIT estimator output caused the fuel flow integrator to 
switch to fan speed control. At the initiation of the deceleration, the estimator output 
had not risen to the point where control of FTIT would have been restored at the inter- 
mediate thrust condition. At the beginning of the deceleration, the fuel flow integrator 
was again clamped since the positive fan speed error would have increased the amoimt 
of downtrim. At the 11. 4-second mark, integral control of fan speed was restored. 
Figures 37(d) and (e) show the responses of the LQR and integral control outputs of noz- 
zle area, respectively. During the initial portion of the acceleration, the trim integrator 
was allowed to uptrim to match the scheduled value of the AP/P parameter. Once the 
AP/P parameter exceeded the scheduled value, however, the integrator was clamped. 
The integrator was again released at the 1. 4-second mark to provide steady-state 
trimming of the nozzle area. The LQR output of the nozzle area acted to reduce the 
state deviations of the fan speed, compressor speed, and afterburner pressure. 

Figures 37(f) to (i) show the LQR and integral control outputs for the fan inlet guide 
vanes and compressor stator vanes. For both sets of vanes, the integrators were used 
to maintain the steady-state vane positions on their nominal schedules. In both cases, 
less than 1. 0° of trim was required. The CIVV integrator was clamped at the 30 per- 
cent thrust condition due to the vanes being fully cambered. Similarly, the RCW in- 
tegrator was clamped at the intermediate thrust condition due to the vanes being fully 
axial. The LQR output of CIW was used to reduce the state deviations of fan speed 
and afterburner pressure. Nearly 6. 0° of downtrim (toward cambered) was used during 
the acceleration and 4. 0° of uptrim (toward axial) during the deceleration. The LQR 
output of RCVV was used to reduce the state deviations of fan speed, compressor 
speed, and afterburner pressure. In this case, nearly 3. 0° of downtrim were required 
during both the acceleration and deceleration. 

Figures 37(j) and (k) show the LQR and integral control outputs of the bleed flow. 

The bleed flow integrator was used to close the bleed valve in steady state. For the 
entire acceleration and steady- state running at intermediate thrust, the negative LQR 
output (positive corresponds to outflow) caused the bleed flow integrator to be clamped. 
Only after the deceleration had reached steady state at the 11. 4-second mark was the 
integrator released to balance the LQR output. The LQR output was used to reduce the 
state deviations of burner pressure, fan speed, and compressor speed. Figures 37(1) 
to (p) show the resulting state deviations. The fuel flow deviation shown in figure 37(p) 
is indicative of the slow increase to the final thrust that was shown in figure 36(f). 

The preceding transient test demonstrated the ability of the multivariable to accel- 
erate the FlOO engine simulation to intermediate thrust while preventii^ a turbine over- 
temperature. The results shown in figures 36 and 37 indicated that the associated con- 
trol logic performed as intended. However, another limit which must not be violated is 
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the minimum burner pressure limit. This limit is a factor at high altitude, low Mach 
number conditions. To demonstrate the multivariable control's ability to maintain suit- 
able levels of burner pressure, the simulation was rim at the 13. 72 km/0. 9 condition. 
The PLA was ramped from 30° to 83° and back to 30° after 10 seconds. At the 30° 
setting, it had been demonstrated that the control uptrimmed the PLAM to 38. 1° to 
maintain the burner pressure at 33. 09 newtons per square centimeter (see fig. 32(c)). 
The transient test was intended to show that the limit would not be violated during the 
rapid deceleration to idle power. Figure 38 shows the results of that test. Figure 38(b) 
shows that the acceleration schedule of fuel flow was, again, not a factor. Also, a slight 
dip in the fuel flow response was again observed. The nozzle area opened up to about 
0. 32 square meter durii^ the acceleration as shown in figure 38(c). At this flight con- 
dition, both the CIVV and RCW went fully axial at the intermediate thrust condition 
as shown in figures 38(d) and (e), respectively. It should be noted that the RCW 
schedule was extended to permit actuator commands of 6. 0° to overcome vane linkage 
hysteresis as indicated in table VI. Figure 38(f) shows that about 0. 5 kilogram per sec- 
ond of bleed flow was utilized during the deceleration. This represented about 3 percent 
of the compressor inlet airflow. Figures 38(g) and (h) show a slight fan speed overshoot 
during the acceleration and a slight compressor speed overshoot during the decelera- 
tion, respectively. The burner pressure limit was maintained during the deceleration 
as shown in figure 38(i). A minimxim burner pressiire of 32. 4 newtons per square centi- 
meter was observed. The dead-band on the burner pressure error driving the fuel flow 
integrator was ±0. 69 newton per square centimeter. Figures 38(b) and (i) show that the 
fuel flow was rapidly uptrimmed once the minimum burner pressure limit was reached. 
Figure 38(j) shows that the FTIT limit was not reached during this test. The estimator 
output matched the actual temperature quite well at the idle condition and during the ac- 
celeration. However, the estimator output was observed to differ from the actual tem- 
perature after the deceleration. Similar behavior was observed for the 3. 048 km/0. 9 
transient as shown in figure 36(j). This was not considered a serious problem, how- 
ever, since the estimator did not perform a control fxmction at the low power settings. 

The corresponding thrust response is shown in figure 38(k). Ninety percent of the 
thrust increase was achieved in 2. 0 seconds with a slow increase to the final thrust 
caused by the dip in fuel flow at the 3. 0-second mark. A thrust overshoot was observed 
durii^ the deceleration and was attributed to the fuel flow uptrim at the low burner 
pressure condition. 

No loss of fan surge margin was observed durii^ the acceleration as shown in 
figure 38(i). An acceptable value of 0. 155 was observed during the deceleration. Dur- 
ing the acceleration, the compressor surge margin dropped to 0. 055 as shown in fig- 
ure 38(m). This was slightly lower than the 0. 075 value observed during hybrid testing 
of the BOM logic but still larger than the 0. 05 minimum requirement. 

Transient tests were also run at sea- level/static, sea-level/1. 2, and 3. 048/1. 2 
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flight conditions with commanded thrust changes larger than 60 percent. In general, the 
multivariable control produced transient behavior similar to that shown in figures 36 
to 38. Table vn summarizes the results of the gross transient tests. The multivariable 
control satisfied the sea-level/static thrust response requirement and produced com- 
parable response times at the other flight conditions. 

In addition to the gross transients, a nximber of other power lever transients were 
run to ensure that the multivariable control logic and implementation worked properly. 
These tests included PLA "bodie” movements, medimn and small amplitude perturba- 
tions of the PLA, and cyclic PLA movements designed to confuse the control logic. 

Figure 39 shows the results of a power lever bodie between the intermediate and’ 

40 percent thrust settings at the 13. 72 km/0. 9 condition. Figure 39 shows that the con- 
trolled engine response tracked the PLA movement quite well. No appreciable speed 
overshoots were observed as shown in figures 39(g) and (h). The FTIT estimator out- 
put, shown in figure 39(i), was observed to peak and remain below the actual tempera- 
ture. This prevented downtrim of fuel flow during acceleration and resulted in the ob- 
served overshoot in FTIT. The FTIT limit was not reached, however, and the over- 
shoot was not considered a serious problem. The thrust responded well to the PLA 
movement with no overshoots as shown in figure 39(j). There appeared to be about a 
0. 4-second lag in the response. Figures 39(k) and {L) show that the fan and compressor 
surge margins did not fall below the minimum requirements for this transient. PLA 
bodies were also run at other flight conditions as indicated in table V. The results of 
those tests were judged to be satisfactory. In general, the engine simulation responded 
well and tracked the PLA movement accurately with no overshoots in speeds, temper- 
atures, thrust, etc. , at those conditions. 

Various medium amplitude movements of the PLA were run to test the LQR and 
integral control gain scheduling algorithm at selected flight conditions. In general, the 
gain switchii^ was smooth and the resultii^ responses were judged satisfactory. How- 
ever, an oscillatory response was observed for an acceleration from PLA = 50° to 
PLA = 70° at the 3. 048 km/0. 9 condition. Figure 40 shows the results from that tran- 
sient. Figures 40(b) and (c) show the responses of fuel flow and compressor speed, re- 
spectively. It appeared that the LQR gains relating the compressor speed deviation to 
fuel flow may have been too high at this condition, which was not one of the LQR design 
points. This problem should be investigated and eliminated prior to engine testing of the 
multivariable control. F^ure 40(d) shows the corresponding response of thrust. Ninety 
percent of the thrust change was achieved in 0. 9 second. Both the fan and compressor 
surge margins were maintained at safe levels for this transient. 

To test the small perturbation response characteristics of the multivariable control, 
±3° PLA snaps were run at several flight conditions as indicated in table V. Prior to 
making the control modifications listed in table VI, data were obtained at sea-level/ 
static, 3, 048 km/0. 6, and 9. 144 km/0. 9 flight conditions. Those results indicated that 
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the small perturbation response characteristics were, in general, satisfactory. Re- 
sponse times ranged from 0. 5 second (at high altitude/high power conditions) to 2. 0 sec- 
onds (at sea-level/static- idle power). However, an underdamped thrust response char- 
acteristic was observed in the 52*^ to 55° PLA range for all three flight conditions. The 
gain scheduling breakpoints were subsequently modified. Small perturbation tests were 
then conducted at the 3. 048 km/0. 9, 13. 72 km/0. 9, sea-level/1. 2, and 3. 048 km/1. 2 
conditions. Results from those tests indicated that the underdamped thrust response 
characteristics had been corrected. Response times ranged from 0. 4 to 1. 0 second for 
these conditions. These response times were judged acceptable relative to the small 
perturbation response requirement of 1. 2 seconds. The small perturbation response 
times are svunmarized in table Vn. 

Figure 41 shows results from the 80° to 83° PLA snap at the 3. 048 km/0. 9 condi- 
tion. A slight fan speed overshoot was observed as shown in figure 41(b). The FTIT 
estimator peak provided fuel flow downtrim and limited the FTIT as shown in figure 
41(d). The 9° overtemperature was ju(%ed acceptable relative to the allowable 27.8° 
overtemperature for 0. 5 second. The thrust response is shown in figure 41(e). Ninety 
percent of the thrust increase was achieved in 0. 5 second. 

Figure 42 shows similar results for the 3. 048 km/1. 2 condition. Figures 42(b) 
and (c) show about 30 rpm overshoots in both rotor speeds. The FTIT estimator peak 
exceeded the gas temperature by 27. 8 kelvins and the limit by 16. 7 kelvins as shown in 
figure 42(d). This caused a premature downtrim of fuel flow and a resulting dip in the 
thrust response as shown in figure 42(e). Ninety percent of the final thrust change was 
achieved in 0. 6 second, prior to the dip in thrust. 

Cyclic power lever movements were established by Pratt & Whitney (ref. 9) and 
used in their preliminary evaluation of the multivariable control logic. This type of test 
was run to ensure that the control logic could not be fooled by cyclic operation around a 
gain switching point. Three such maneuvers were selected for the hybrid evaluation and 
programmed on an analog function generators. The cyclic PLA tests were conducted 
at the flight conditions listed in table V. Figure 43 shows a typical result from that in- 
vestigation. The PLA movement shown in figure 43(a) was run at the 13. 72 km/0. 9 
condition. Figures 43(b) to (f) show the resulting engine responses. Figure 43(d) shows 
that the thrust tracked the PLA movement quite well with little phase lag. The fan and 
compressor surge margins were maintained at safe levels as shown in figures 43(e) 
and (f), respectively. Similar results were obtained at the other flight conditions. 


Control Flexibility - Fast Acceleration Design 

The results of Pratt & Whitney’s preliminary evaluation of the multivariable control 
logic (ref. 9) indicated that the potential existed for achieving faster thrust response. 
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The previously discussed results from the hybrid evaluation also indicated this fact. 

While fast thrust response was not a specific objective of the MVCS program, it was 
felt that achieving faster response through a straightforward design iteration without 
changing the structure of the control wovild demonstrate the flexibility of the design ap- 
proach and resulting implementation. 

The Systems Control Inc. approach to achieving faster response with the multivari- 
able control was to redesign the transition control rate limits. These rate limits define 
the transient model which the regulator attempts to follow. For the hybrid evaluation, 
the redesign was accomplished at sea-level/static conditions. The resulting rates were 
implemented in the SEL810B computer and were evaluated for PLA snaps from the idle 
to the intermediate thrust setting and from the 30 percent to the intermediate thrust 
setting. The resulting transient data were compared with the corresponding results ob- 
tained with the normal rate limits. Figure 44 shows the results of that comparison for 
the 30 percent to intermediate thrust transient. Figure 44(b) shows that the fast- 
acceleration rates resulted in the fuel flow approaching the acceleration limit (fig. 17) 
during the acceleration. The resulting responses of fan speed, compressor speed, FTIT, 
and thrust are shown in figures 44(c) to (f), respectively. No overshoots were observed 
during the acceleration and only slightly more overshoot was observed during the decel- 
eration with the fast-acceleration rates. The FTIT estimator output caused fuel flow 
downtrim at the 0. 6-second mark resulting in the dip in the FTIT response shown in 
figure 44(c). The faster thrust response, shown in figure 44(f), exhibited a 1. 4 -second 
response time as compared with the normal response time of 2. 2 seconds. The faster 
thrust response was achieved with little or no loss in fan or compressor surge margin 
as shown in figures 44(g) and (h), respectively. 

Similar results were obtained for the PLA snap from the idle to the intermediate 
thrust setting. A reduction in the response time from 3. 2 to 2. 2 seconds was achieved 
with a loss of 0. 012 in fan surge margin and a loss of 0. 07 in compressor surge margin. 


Afterburner Ignition Tolerance 

Although LQR control of afterburner fuel flow was not a part of the multivariable 
control task, the control requirements dictated that the control be tolerant of external 
disturbances such as afterburner ignition pulses. Figure 10 specified the maximum al- 
lowable amplitudes of afterburner pressure pulses. These values were based on BOM 
control results for which the exhaust nozzle is opened in anticipation of the afterburner 
ignition. It was hoped that the LQR would provide sufficient regulation based on the 
sensed afterburner pressure deviations. 

Afterburner ignitions were simulated by a step-wise increase in the afterburner fuel 
flow. The amount of injected fuel was based on the current afterburner schedule and the 
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intermediate value of burner pressure for the fl^ht condition being studied. Figure 45 
shows the engine and control variable responses at the 13. 72 km/0. 9 condition. Fig- 
ures 45(b) to (d) show the resultant deviations in fan speed, compressor speed, and after- 
burner pressure, respectively. The amplitude of the afterburner pressure spike was 
approximately 7. 2 percent of the intermediate pressure. This was about 2 percent 
higher than the desired amplitude. Figures 45(e) and (f) show that losses of 0. 065 in fan 
surge margin and 0. 035 in compressor surge margin were observed. Figures 45(b) 
and (g) show the result of the integral control action which attempts to maintain the fan 
speed and the fan discharge AP/P parameter at their scheduled values. While one 
would expect the suppression of the pressure spikes to be accomplished by opening the 
nozzle area, the LQR, in fact, opened the nozzle only slightly (see fig. 45(h)) while mov- 
ing the fan inlet guide vanes 2® cambered. F^ure 46 shows the results of the after- 
burner ignition tests at all of the evaluated flight conditions. These results indicated 
that the LQR design, without anticipatory nozzle opening, did provide pressure spike 
suppression. It appears that incorporating an extended nozzle area schedule and the 
necessary prefill logic for the afterburner fuel manifolds would be compatible with the 
multivariable control logic. The LQR action would, in that case, supplement the nozzle 
schedule in suppressir^ the pressure spikes. 


Sensor Failure Accommodation 

Althoxigh sensor failure accommodation was not a requirement for the multivariable 
control design, hybrid tests were conducted to identify critical sensors in order that ap- 
propriate software provisions could be made prior to full-scale engine tests. Various 
modes of sensor failure were studied at the 9. 144 km/0. 9 condition and an intermediate 
thrust setting. Both full-scale (saturation) and sensor loss tests were conducted. The 
results of the failure study are summarized in table Vm. The control structure with 
its engine protection logic proved adequate for most failures. Two critical failures 
were identified, however. They were saturation of the fan inlet total pressure (P 2 ) sen- 
sor and the loss of the fan speed sensor. Figure 47 shows the effects of the Pg sensor 
saturation. The high P 2 signal resulted in a sudden increase in the scheduled values 
of fuel flow, burner pressure, and afterburner pressure. The feed-forward path for 
the scheduled control values drove the fuel flow command up as shown in figure 47(b). 
The 680 kilograms per hour downtrim limit on the fuel flow integrator prevented it from 
limiting the fan speed which exceeded its limit at about the 9. 5-second mark as shown in 
figure 47(c). Similarly, the compressor speed exceeded its limit at about the 8. 5- 
second mark as shown in figure 47(d). The FTIT estimator is dependent on having an 
accurate schedule of FTIT. The FTIT schedule is, however, a function of the sensed 
fan inlet air density. The high P 2 signal caused an erroneous estimator output as 
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shown in figure 47(e). In spite of the overtemperature at the 1. 4-second mark, no down- 
trim of fuel flow was requested. Facility shutdown logic, based on the sensed FTIT, 
would not be active until 4. 5 seconds after the sensor failure due to the slow sensor dy- 
namics, thus pointing out the need for some form of Pg sensor failure detection logic. 
Figures 47(f) through (h) show the correspondii^ thrust and surge margin responses. 

The loss of fan surge margin also indicates the necessity of some form on detection and 
accommodation logic for the sensor. The ?2 sensor saturation test was subse- 
quently repeated with the downtrim limit removed from the fuel flow integrator. Nearly 
2700 kilograms per hour of downtrim was requested but could not prevent the overtem- 
perature and fan overspeed from occurring. The compressor speed was maintained 
within limits, however. 

The loss of the fan speed sensor resulted in increased fuel flow and nozzle area 
trims from the regulator. The fuel flow integrator was also driven to its uptrim limit 
by the false fan speed deviation. This resulted in an overtemperature at the 0. 2-second 
mark, a fan overspeed at the 0. 28-second mark, and a compressor overspeed at the 
0. 57-second mark. Based on these results, it is clear that detection and accommodation 
Ic^ic are required for the fan speed sensor also. 


Flight Condition Transition 

The results from the previously described tests indicated that the multivariable con- 
trol would provide satisfactory regulation at each of the evaluated flight conditions. 
However, it was noted that the FTIT estimator could be confused by sensor noise and/ 
or changes in the flight condition. In changing flight conditions during the hybrid evalua- 
tion, no attempt was made to control the rates of change of the fan inlet conditions (P 2 
and T 2 ). Therefore, a test was conducted which simulated representative aircraft ac- 
celerations and decelerations as defined in the CONTROL REQUIREMENTS section. At 
sea level, the flight Mach number was ramped from 0. 0 to 0. 6 in 32 seconds as shown 
in figure 48(a). This resulted in the fan inlet pressure and temperature responses shown 
in figures 48(b) and (c), respectively. The peak rates of change were 0. 17 newton per 
square centimeter per second and 1. 3 kelvins per second. Figures 48(d) to (f) show the 
resulting responses of fuel flow, FTIT, and thrust. Figure 48(d) shows that the fuel flow 
was not downtrimmed adequately until late in the acceleration. Although the fan and 
compressor speeds did not exceed their limits, the FTIT did exceed its limit by as 
much as 11. 1 kelvins as shown in figure 48(e). The FTIT overtemperature indicated 
that the estimator output did not respond accurately to the FTIT sensor output and the 
commanded fuel flow. A slight thrust overshoot was also observed as shown in fig- 
ure 48(f). 

A corresponding sea-level deceleration was simulated with the Mach number ramped 
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from 0. 6 to 0. 0 in 11 seconds as shown in figure 49(a). This resulted in maximum 
rates of change of PT2 and TT2 of -0. 59 newton per square centimeter per second 
and -2. 8 kelvins per second as shown in figures 49(b) and (c). For this deceleration, 
the estimator provided adequate downtrimming of the fuel flow to maintain FTIT within 
limits as shown in figures 49(d) and (e). The results of the flight condition transition 
test point out the need to take a closer look at the estimator behavior prior to the engine 
tests. 


SUMMARY OF RESULTS 

The primary objective of the evaluation was to verify the multivariable control 
logic and its implementation so as to ensure safe and stable operation of the FlOO en- 
gine during subsequent altitude tests. The results of the evaluation indicated that the 
control logic and its implementation on the SEL810B digital computer will be capable of 
controlling the engine throughout its operating range. The specified engine limits were 
not violated during normal steady- state and transient operation. 

The results of the steady-state evaluation showed that the reference point schedules 
matched the simulation model of the engine at most operating conditions. Where dis- 
crepancies were noted, they did not affect the normal control action. An observed dis- 
crepancy in the FTIT schedule at the 19. 81 km/0. 9 condition, coupled with the low 
burner pressure at that condition, did result in an overtemperature whenever changes 
in the fan inlet conditions caused resetting of the FTIT estimator to the low scheduled 
value. A subsequent change in the FTIT estimator logic to prevent the reset action at 
the low burner pressure condition eliminated this problem. 

The hybrid evaluation indicated that the multivariable control could accommodate 
the slight differences between the hybrid model and the CCD1103-1. 0 simulation. While 
the control should be able to handle expected engine-to-engine differences, it is recom- 
mended that the reference point schedules be trimmed as soon as FSER engine data 
becomes available. 

Simulation errors, caused by simplified thrust and surge margin calculations, did 
make it difficult to quantitatively evaluate the steady-state, closed-loop performance of 
the controls. However, previous studies involving a digital implementation of the BOM 
control Ic^ic and the hybrid simulation provided an opportunity to ’•adjust’* the MVC 
data. The adjusted values of thrust indicated that the multivariable control matched 
the CCD1103-1. 0 intermediate thrust values (within the specified 1 percent) at all but 
the siQ)ersonic flight conditions. This result is substantiated by data provided by SCI 
from their version of the CCD1103-1. 0 simulation with FORTRAN coding of the multi- 
variable control logic. Those data also showed low thrust values at the si 5 )ersonic flight 
conditions. This is attributed to the fact that the supersonic flight points were not fac- 
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tored into the reference point schedules. It is recommended that adjustments (if pos- 
sible) be made to the reference point schedules to improve the h^h Mach number per- 
formance. 

Authority lim its on the fuel flow and nozzle area trim integrators were implemented 
to provide protection in the event of sensor failures. It was determined, however, that 
the assumed limits of ±0. 028 square meter and ±680. 4 kilograms per hour did prevent the 
steady-state matching of the schedules at three supersonic flight conditions. This was 
caused by the aforementioned scheduling errors at the supersonic flight conditions. 

While no engine limits were violated due to the trim saturation, the schedules should be 
refined prior to the engine tests. 

The proportional (LQR) plus integral control structure provided good fan operating 
point control over the range of ei^ine operation. The integral control provided adequate 
limiting of fan speed, FTIT, and burner pressure at all flight conditions. 

The LQR and transition control resulted in satisfactory transient responses at most 
operating conditions. The gain switchii^ for gross transients was smooth. The break 
between low and high power gains was shifted to a higher power condition to eliminate 
underdamped behavior at low power settings. Unsatisfactory response characteristics 
were noted at only two conditions. Oscillatory behavior was observed during a 50° to 
70° PLA snap at the 3. 048 km/0. 9 condition. This was attributed to a high regulator 
gain relatii^ fuel flow to the compressor speed deviation at this condition. As described 
in the MULTIVARIABLE CONTROL section, the LQR gains at this condition are de- 
pendent on the gain scheduling algorithm since none of the stored LQR gain matrices 
were designed at this flight condition. It is recommended that this behavior be cor- 
rected prior to the engine tests. Also, an excessive FTIT estimator overshoot was ob- 
served during an 80° to 83° PLA snap at the 3. 048 km/1. 2 condition. This resulted in 
a premature downtrimming of fuel flow during the acceleration and a corresponding dip 
in the thrust response. FTIT anomalies were observed durii^ various phases of the 
evaluation. In addition to the aforementioned overshoot, a slowly diminishii^ error be- 
tween the estimator output and the actual gas temperature was observed at the conclu- 
sion of many of the transients. Although the estimator fulfilled its basic function of 
providing a form of anticipation of impending overtemperatures, it is recommended that 
the estimator design be refined prior to the engine tests so as to provide predictable 
behavior durii^ transients and flight condition transitions. 

The specified response requirements were satisfied for all small and large ampli- 
tude transients with the exception of the small (3°) PLA snaps at the sea-level/static, 
idle condition. It should be noted that the 1. 2-second response time requirement for 
small perturbations was adopted due to a lack of specificity in the requirements. 

The flexibility of the control structure and the design methods was demonstrated by 
implementing a set of fast-acceleration rate limits in the transition control. The re- 
sulting sea-level/static acceleration from idle to intermediate thrust showed a reduction 
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in the response time from 3. 2 to 2. 2 seconds. This compafed with the current BOM- 
controlled response time of 2. 8 seconds. 

Afterburner ignitions were simulated to test the multivariable control's tolerance 
to external disturbances. The results of those tests indicated that the LQR design did 
provide pressure spike suppression at all flight conditions. The LQR was observed to 
use both the exhaust nozzle area and the fan inlet guide vanes (CIVV) to suppress the 
pressure spikes. 

The results of the sensor failure study of the 9. 144 km/0. 9 condition indicated that 
most sensor failures would result in a safe, downtrimmii^ to a part -power condition. 
The saturation of the Pg sensor or the loss of the fan speed sensor, however, resulted 
in a catastrophic overspeed and overtemperature condition. Therefore, it is recom- 
mended that some form of sensor failure detection and accommodation logic be imple- 
mented in the multivariable control if redundant sensors cannot be provided for the en- 
gine tests. It is further recommended that some form of the current acceleration 
schedule of 4'^^4)m implemented since it provides fuel flow cutback in the 
event of a loss of the compressor speed signal or a compressor overspeed. The results 
of the hybrid evaluation showed that the limiting of the reference point schedule of com- 
pressor speed is not sufficient to prevent overspeeds from occurring althov^h normal 
engine operation would seem to preclude this from happening. 

The results of the simulated aircraft accelerations and decelerations at sea level 
showed that the FTIT estimator Ic^ic could be confused by changing fan inlet conditions. 
This resulted in overtemperatures in some cases. Those results point out the need to 
take a closer look at the estimator design prior to the engine tests. 


CONCLUDING REMARKS 

The results of the evaluation indicate that the engine testii^ of the multivariable 
control should be conducted as planned. The multivariable control produced satisfactory 
performance for all but a few test conditions. With minor modifications, the control 
should provide acceptable steady-state and transient performance at all flight conditions 
and power settings. 

Considering the fact that in approximately 10 months a control design for the FlOO 
engine was accomplished, it must be concluded that the computer-aided approach to 
designing multivariable controls is a practical solution to the ermine control problem. 

The results of the evaluation indicate that the Systems Control Inc. approach to solvii^ 
the nonlinear control problem is also practical. The implementation of the control logic 
on the SEL810B digital computer has demonstrated that a practical LQR control can be 
implemented on a digital computer having many of the characteristics of flight-qualified 
computers. The core requirements (7. 1 K) and update time (10 msec) are comparable 
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to the requirements of flight-qualified digital controls such as the IPCS control and the 
FlOO electronic supervisory control. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 22, 1977, 

505-05. 
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TABLE I. - FlOO TRANSIENT THRUST 
RESPONSE TIME REQUIREMENTS 


[All times for sea-level, static, standard day with no bleed 
flow or power extraction.] 


Thrust 

change, 

percent 

Power lever angle change 

Idle to in- 

30 Percent to 

n 

Intermediate 

Idle to 

termediate 

intermediate 

to idle 

maximum 



Time, 

sec 


90 

4.0 

3.5 

3.0 

8.0 

98 

15.0 

15.0 

20.0 

15.0 
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TABLE n. - MULTIVARIABLE CONTROL 
COMPUTER SPECIFICATIONS 


SEL810B digital computer 

Two 16 -bit accumulators 
Memory specifications: 

24K magnetic core 
Cycle time, 0. 75 jusec 
Expandable to 32K 

Two’s -compliment, fixed-point multiply and divide: 
Add time, 1. 5 jLisec 
Multiply time, 4. 5 iisec 
Divide time, 8. 25 jusec 
Double precision arithmetic 
Infinite indirecting 
Infinite indexing 
Direct memory access 
28 Levels of vectored priority interrupt 
66 Total instructions 

Analog acquisition unit®' 

Two multiplexors sample and hold circuits and ana- 
log to digital converters 
64 Input channels for each multiplexor 
Input voltage range, ±10 V 

12 Bit + sign data resolution (two's compliment) 
Digitizing rate, 50 ju sec/sample 
Percent error with calibration, 0. 073 

Analog output unit®^ 

26 Digital to analog conversion channels: 

10 to 12 bit (plus sign) 

16 to 11 bit (plus sign) 

Output voltage range, ±10 V 
Slew rate, 1 V/sec 

Tektronix 4010A peripheral system®^ 

Tektronix 4010 scope terminal: 

1600 Baud capacity 
RS232 and teletype interface 
Paper tape reader punch 
Floppy disk: 

Store 262, 144 bytes/disk - ^^660 meters of paper 
tape 

Write speed, 400 bytes/sec 


®SEL810B peripherals. 


TABLE m. - CORE REQUIREMENTS FOR 
MULTIVARIABLE CONTROL 


(MVC) PROGRAM 


Software 

Core 

storage, 

words 

Block data: 


Schedules 

1 226 

Matrices 

450 

Total 

1 676 

MVC program: 


Main 

1 305 

FXrr estimator 

315 

Set point 

547 

Gain control 

830 

Transition control 

634 

Integral control 

776 

LQR control 

332 

Engine protection 

202 

Function generation 

518 

Total 

5 459 

General purpose input/output and debug routines: 


INFORM and DATAO and SAMPLE 

3 613 

Input/output drivers 

674 

Floating point software 

519 

Total 

4 806 

Grand total 

11 941 




TABLE IV. - FIDO RESEARCH 


SENSOR CHARACTERISTICS 


Sensor 

Characteristics 

Fan speed 

1 

(0. 03 S + 1) 

Compressor speed 

1 

(0.05 S+ 1) 

Pressures 

1 

(0. 05 S + 1) 

Faa inlet temperature 

1 

(1. 5 S + 1) 

Fan turbine inlet tempera- 
ture 

0.31 1 0.69 

(0.6 S+l) (5.5 S+1) 

Fan discharge temperature 

(3S + 1)(15. 5 S + 1) 

(tS + 1)(20S + 1)(0. 6S + 1) 


where 


T = f(Wj3) 


TABLE V. - FIDO MULTIVARIABLE CONTROL SYNTHESIS PROGRAM- 


HYBRID EVALUATION TESTS 


Altitude (km) 

Steady state 

Power lever angle 

Disturbances 

Sensor 

Mach number 

PLA = 83° 

PLA<83° 

(PLA) transients 

After- 

Maneuver 

failures 


Large 

Small 

Cyclic 

burner 


0/0 

a 

• 

• 

• 

• 

• 

• 

• 


3.048/0.6 

• 

• 

• 

• 

• 

• 

b 

o 


3. 048/0. 9 

• 

• 

• 

• 

• 

• 

o 


9. 144/0. 9 

• 

• 

• 

• 

• 

• 

o 

• 

13. 72/0. 9 

• 

• 

• 

• 

• 

• 

o 


19. 81/0. 9 

• 





• 

o 


0/1.2 

• 

• 

• 

• 


• 

o 


3.048/1.2 

• 

• 

• 

• 


• 

o 


6.096/1.8 

• 





• 

o 


22. 86/1. 8 

• 





• 



17. 83/2. 15 
12. 19/2. 2 

• ! 
• 





• 

o 


19. 81/2. 5 

• 





• 

o 



Specific test, •. 

'^Evaluated during flight condition transients between tests, o. 



TABLE VL - PREEVALUATION CONTROL MODIFICATIONS 
Symptom I Remedy 


Under damped behavior at midpower 
for cyclic PLA inputs 

Erratic FTIT estimator behavior 

RCW off- schedule at high-power 
due to actuator hysteresis 

Bleed command noise 

Excessive nozzle area downtrims 
at 19. 81 km/0. 9 condition due to 
sensitivity to A P/P calculation 
accuracy 

Anticipated problems caused by 
sensor failures or scheduling 
errors 


Extension of gain transition region to 
lower gains in midpower range 

Dead- zone placed in estimator initiali- 
zation logic 

Extension of schedule at high compres- 
sor speeds 

Negative bleed command bias with 
added downstream limits 

Authority limit placed on nozzle area 
trim 

Limits placed on fuel flow trim and 
trim rate 


TABLE Vn. - POWER LEVER ANGLE (PLA) 
RESPONSE RESULTS 


[Specified response time requirements in parentheses, A 1. 2- 
sec response time requirement assumed for ±3® PLA 
changes.] 


Altitude (km) 
Mach munber 


PLA change, deg 


20-23 

30-33 

52-55 

80-83 

Idle-83 

35-83 


Thrust response times, sec 


0/0 

2.0 

1.4 

0.8 

0.7 

3. 2 (4. 0) 

2. 2 (3. 5) 

3. 048/0. 6 

1.2 

1.0 

.6 

.6 

— 

— 

3. 048/0. 9 

1.0 

1.0 

.8 

.5 

3.0 

2.1 

9. 144/0. 9 

— 

.7 

.5 

.5 

— 

— 

13. 72/0. 9 

— 

— 

.4 

.4 

2.0 

— - 

0/1.2 

--- 

— 

.6 

.5 

1.5 

— 

3. 048/1. 2 

— 

— 

.75 

.6 

1.5 

— 
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TABLE Vin. - RESULTS OF HARD SENSOR 


FAILURES ON MULTIVARIABLE 
CONTROL AT 9. 144 km/0. 9 
AND PLA = 83° 


Sensor 

Mode 

Result 

MN 

Loss 

None 


Saturation 

Slight downtrim 

^2 

Loss 

Downtrim to 38 percent thrust 


Saturation 

Overtemperature, overspeed^ 

^2 

Loss 

None 


Saturation 

Downtrim to 63 percent thrust 


Loss 

Overtemperature, overspeed^ 


Saturation 

Downtrim to 60 percent thrust 

N2 

Loss 

Downtrim to 34 percent thrust 


Saturation 

Downtrim to 58 percent thrust 

P4 

Loss 

Downtrim to idle 


Saturation 

None 

^6 

Loss 

None 

Saturation 

Downtrim to 74 percent thrust 

(AP/P)i3 

Loss 

Slight downtrim 

Saturation 

None 

'"f,4 

Loss 

None 


Saturation 

Downtrim to 68 percent thrust 

^Catastrophic failure. 




Figure 1. - Multivariable Control Synthesis program - design and evaluation process. Pratt & Whitney 
Aircraft. P&WA; Systems Control inc., SCI. 









Fan corrected airflow, 


Cambered 




Figure 7. - Fan corrected airflow limits. 




Figures. -Typical power extraction requirements: standard 
day. 



Figure 9. -Typical bleed requirements: standard day, all 
power settings. 



Figure 10. - Afterburner ignition pressure pulse 
requirements. 



PLAM 

MN 

h 



Figure 11. - R 00 Multivariable Control System. 



Figure 12, - Reference point scheduling algorithm. 
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Rate- limited 
scheduled 
variables, 
u 

^nom* ^nom 


Figure 13. - Transition control algorithm. 
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state d^iations 





Fan 

speed 

Compressor 

speed 

Afterburner 

pressure 

Fuel 

flow 

Burner 

pressure 

Fuel 

flow 

X 

X 

X 

X 

X 

Nozzle 

area 

X 

X 

X 

0 

0 

CiVV 

X 

0 

X 

0 

0 

RCW 

X 

X 

X 

0 

0 

Bleed 

flow 

X 

X 

X 

0 

X 


Figure 14. - LQR gain matrix structure. 


Controls < 


State deviations Spares Integral output deviations 



f 

Fan 

speed 

Compressor 

speed 

Afterburner 

pressure 

Fuel 

flow 

Burner 

pressure 


r 

AP/P 

parameter 

CIW 

RCW 

Bleed 

flow 

Fan 

speed 

FTIT 

maximum 

limit 

Burner 

pressure 

minimum 

limit 



Burner 

pressure 

maximum 

limit 

Fuel 

flow 

X 

X 

X 

X 

X 


X 

0 

0 

0 

X 

X 

X 

X 

Nozzle 

area 

X 


X 

0 

0 


X 

0 

0 

0 

X 

X 

X 

X 

CIVV 

X 

0 

X 

0 

0 


X 

X 

0 

0 

X 

0 

0 

0 

RCW 

X 

X 

X 

0 

! 0 


0 

0 

X 

0 

0 

0 

0 

0 

Bleed 

flow 

X 

X 

0 

0 

X 


0 

0 

0 

X 

0 

0 

0 

0 


Figure 15. - Gain matrix structure. 



Figure 16. - Gain schedule algorithm. 
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Figure 17. - Engine protection algorithm. 
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Figure 18. - Fan turbine inlet temperature estimator. 
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Figure 25. - Schematic representation of exhaust nozzle actuator system. 
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Figure 26. - Schematic representation of CIW actuator system. 
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Figure 27. - Schematic representation of fuel valve actuator system. 



Figure 28. - Schematic representation of RCVV actuator system. 
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Scheduled burner pressure, N/cm^ 
(c) Burner pressure. 



Scheduled afterburner pressure, N/cm^ 
(d) Afterburner pressure. 


Figure 29. - Steady-state, uninstalled performance of FlOO real-time simulation with multivariable control. Sea-level, static, standard 
day conditions. 
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(g) Adjusted net thrust. 





Scheduled fan discharge AP/P parameter 


(f) Fan discharge AP/P parameter. 



Power lever angle, deg 


(h) Adjusted fan surge margin. 

Figure 29. - Continued. 
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Power lever angle, deg 


(h) Adjusted fan surge margin. 


Figure 30. - Concluded. 


(I) Adjusted compressor surge margin. 
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Figure 31. - Steady-state, uninstalled performance of FlOO real-time simulation with multivariable control. Altitude, 9. 144 kilometers; Mach num 
ber, 0.9; standard-day conditions. 
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(g) Adjusted net thrust. 



(f) Fan discharge AP/P parameter. 
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(h) Adjusted fan surge margin. 

Figure 31. - Continued. 



Fan speed, rpm 



Power lever angle, deg 


(i) Adjusted compressor surge margin. 
Figure 31. - Concluded. 



Scheduled fan speed, rpm Scheduled compressor speed, rpm 


(a) Fan speed. (b) Compressor speed. 

Figure 32. - Steady-state, uninstalled performance of FlOO real-time simulation with multivariable control. Altitude, 13.72 kilo- 
meters; Mach number, 0.9; standard-day conditions. 
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(h) Adjusted fan surge margin. (i) Adjusted compressor surge margin. 

Figure 32. - Concluded. 
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Scheduled FTIT -maximum FTIT, K Scheduled fan discharge AP/P parameter 

(e) Fan turbine inlet temperature. (f) Fan discharge AP/P parameter. 



Power lever angle, deg 

Ig) Adjusted net thrust. (h) Adjusted fan surge margin. 

Figure 33. - Continued. 



FTir - maximum FTIT, K Burner pressure, HIcnr 
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Figure 34. - Continued. 
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Figure 35. - Steady-state, uninstalled performance of FIDO real-time simulation with multivariable control. Standard-day conditions; intermediate 
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ire 35. - Concluded. 


Nozzle area, Power lever angle, deg 



Time, sec 

(c) Nozzle area. (d) Fan inlet guide vanes (CIVV). 


Figure 36. - Simulated FlOO response to 35° to 83° power lever angle (PLA) snap at 3. 048 kilometers and Mach 0. 9. Multivariable control; PLA = 35° 
provides 30'percent thrust. 







Time, sec 

(k) Afterburner pressure. (i) Net thrust. 

Figure 36. - Continued. 



(m) Fan surge margin. 



(n) Compressor surge margin. 



(o) Fan discharge AP/P parameter. 


Figure 36. - Concluded. 

















Fan surge margin FTIT - maximum FTH. 






(c) Compressor speed. 

Figure 40. - Simulated FlOO response to 5(P to 70° power 


Net thrust, percent of intermediate Fuel flow, kg/hr 




angle (PLA) snap at 3.048 kilometers and Mach 0.9. Multivariable control. 









Time, sec 


(e) Fan surge margin. (f) Compressor surge margin. 

Figure 43. - Concluded. 



Time, sec 

(a) Power lever angle. (b) Fuel flow. 


Figure 44 - Comparison of normal and fast FlOO simulation responses to 35° to 83° power lever angle (PLA) snap at sea-level, static conditions. Multi 
variable control. 









Regulator output, Compressor surge margin 












1. Report No. 


2. Government Accession No. 


NASA TP-1056 | 

4. Title and Subtitle pioQ MULTIVARIABLE CONTROL SYNTHESIS 


PROGRAM - EVALUATION OF A MULTIVARIABLE CONTROL 
USING A REAL-TIME ENGINE SIMULATION 

7. Author(s) 

John R. Szuch, James F. Soeder, Kurt Seldner, 
and David S. Cwyn^ 

9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

15. Supplementary Notes 


3. Recipient's Catalog No. 


5. Report Date 

October 1977 

6. Performing Organization Code 

8. Performing Organization Report No. 

E-9170 

_ 10. Work Unit No. 

505^05 

1 1 . Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Paper 

14. Sponsoring Agency Code 


16. Abstract 

The FlOO Multivariable Control Synthesis program is a cooperative effort by the Air Force 
Aero- Propulsion Laboratory and the NASA Lewis Research Center aimed at accomplishing 
the design, evaluation, and testing of a practical, multivariable, linear quadratic regulator 
control for the FlOO turbofan engine. This report covers the NASA evaluation of the multi- 
variable control logic and implementation. The evaluation utilized a real-time, hybrid 
computer simulation of the engine. The results of the evaluation are presented, and recom- 
mendations concerning future engine testing of the control are made. The results indicated 
that the engine testing of the control should be conducted as planned. 


17. Key Words {Suggested by Author(s}) 


18. Distribution Statement 


Simulation; Hybrid computer; Real-time; Unclassified - unlimited 

Turbofan; Multivariable control STAR Category 07 


19, Security Ciassif, (of this report) 

20. Security Ciassif. (of this page) 

1 21. No. of Pages 

22. 

Unclassified 

Unclassified 

103 



For sale by the National Technical information Service, Springfield, Virginia 22161 


NASA-Langley, 1977 



National Aeronautics 
Space Administration 

Washington, D.C. 
20546 

Official Business 
Penalty for Private Use, i 




and 


THIRD-CLASS BULK RATE 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



2 1 10, A, 092377 S00903DS 

DEPT OF THE AIB FCBCE 
AF WEAPONS LABOBAIOSY 
ATTN: TECHNICAL LIBRABY (SUL) 
KIETLAND AFB NM 87117 



POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 




